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1.     INTRODUCTION 

The  purpose  of  this  paper  is  to  explore  the  use  of  lunar  surface 
photography  in  order  to  achieve  the  photogrammetric  transfer  of  available 
selenographic  coordinates  from  future  lunar  landing  sites  to  neighboring, 
photoidentifiable  features.    It  can  be  implied  from  the  procedures  developed 
that  overhead  photography,  were  it  available,  could  be  utilized  and  would 
provide  a  material  strengthening  of  the  total  solution.     By  the  methodic 
selection  of  features  and  confirmation  that  they  can  in  reality  be  identified 
from  orbital  photography,  a  modest  selenodetic  control  system  can  be 
expanded  into  a  net  that  could  ultimately  control  all  future,  manned  or 
unmanned,  orbital  photographic  missions. 


2.     ABSTRACT 

For  centuries  man  has  scrutinized  the  moon  in  one  manner  or 
another  and  postulated  theories  concerning  its  size,   shape,   origin,   and 
other  general  characteristics.    With  the  passage  of  time  and  the  improve- 
ment of  equipment  and  observation  techniques  the  desire  for  more  explicit 
information  concerning  earth's  nearest  celestial  neighbor  has  become  acute. 
In  fact,   as  the  moment  approached  when  man  would  actually  set  foot  on  the 
lunar  surface,  the  need  for  such  information  became  vital.     The  following 
historical  review  briefly  outlines  man's  effort  to  improve  his  knowledge 
in  one  of  the  pertinent  regions  of  selenodesy  —  selenodetic  control. 

The  remainder  of  this  paper  explores  a  method  of  improving  the 
existing  selenodetic  control  by  employing  available  lunar  surface  photography 
supplemented  by  that  obtained  from  lunar  orbit.     Following  the  results  of 
this  experiment  an  ideal  model  is  submitted.     The  unknowns  associated  with 
this  model  are  perturbed  within  realistic  limits  by  a  random  number  gener- 
ation program.    This  provides  a  theoretical  indication  of  the  accuracy  that 
could  be  anticipated  assuming  there  is  reasonable  adherence  to  the  suggested 
procedures. 

Finally,  conclusions  are  drawn  and  reasonable  recommendations  are 
offered  to  improve  selenodetic  control  by  the  photogrammetric  transfer  of 
known  or  assumed, local  or  astronomic  coordinates  of  a  lunar  landing  site  to 
neighboring  features  that  may  be  photoidentified  from  orbital  photographs. 


3.  HISTORICAL  REVIEW 

For  decades  in  the  past  to  the  present  day  the  task  of  surveying  the 
moon  has  engaged  the  efforts  of  many  astronomers.     In  early  1959  the 
launching  of  LUNIK  I  by  Russia,   and  subsequently,    POINEER  IV  by  the 
United  States:    ".  .  .  opened  the  first  modern,  post  telescopic  phase  of  lunar 
exploration.  "  or,   at  least,   introduced  a  tantalizing  new  dimension  [30]. 

During  the  some  seventy  years  prior  to  the  launching  of  the  first 
lunar  space  probes  the  establishment  of  selenodetic  control  was  founded  on 
direct  astronomic  angular  observations  and  indirect  angular  observations 
through  astronomic  photography.     Essentially,   it  was  based  on  heliometric 
observations  which  consist  of  measurements  of  position  angles  and  angular 
distances  between  a  reference  point  on  the  lunar  surface  (Mbsting  A  is  the 
fundamental  point)  and  the  lunar  limb.     Observations  at  mean  libration 
permit  a  best-fit  circle  of  the  lunar  disc  to  be  established.     The  center  of 
this  circle  is  defined  as  the  projection  of  the  origin  of  the  coordinate 
system  (the  dynamic  or  mass  center  of  the  moon)  upon  the  lunar  surface 
and  its  radius  to  be  the  mean  radius  of  the  moon.     Thus,  the  center  of 
figure  is  equated  to  the  center  of  mass  and  in  the  adjustment  of  the  helio- 
metric observations  this  injects  the  so-called  center  of  figure  bias.     The 
adjustment  provides  corrected  values  of  physical  lunar  libration  parameters 
and  the  coordinates  of  the  reference  point  as  well  as  the  mean  radius  of  the 
moon  [26]. 

The  heliometer  was  first  developed  by  Bouger  in  1748  and  later 
modified  by  Dollond.     It  consists  of  a  refractor  telescope  with  two  semi-lenses 
which  may  form  a  single,   superimposed  image  of  two  object  points  at  the 
principle  focus.    The  angular  distance  between  the  two  object  points  formed 
on  the  focus  is  equal  to  the-distance  between  the  centers  of  the  semi-lenses 
when  one  slides  parallel  to  a  line  of  section  upon  the  other  [24].  It  was  used 
to  measure  the  diameter  of  the  moon  at  the  end  of  the  18th  century  byLalande 


and  by  Bessel  in  1839  to  investigate  lunar  physical  librations.     It  was 
Bessel  that  developed  the  procedures  for  measurement  that  remain  basically 
intact  today. 

Heliometric  observations  are  limited  by  the  resolving  power  im- 
posed by  their  relatively  small  aperatures  (4-7  inches).    The  "Rayleigh 

criterion: 

6=1.22% 

9  =  minimum  angle  resolved  in  minutes 

X  =  wavelength  of  light 

D  =  diameter  of  objective  lens 

theoretically  indicates  that  a  six  inch  aperature  provides  a  minimum  resolu- 
tion of  0.  75  arc  seconds  or  well  over  a  kilometer  on  the  moon' s  surface.     A 
further  limitation  is  based  on  atmospheric  refractivity  [20]. 

Selenodetic  control  systems  derived  from  earth-based  lunar  photo- 
graphy generally  rely  heavily  on  heliometric  observations.    The  reduction 
of  these  observations  provide  the  libration  parameters  (f,  I)  and  the  coor- 
dinates of  fundamental  points.     These  reference  points  provide  the  orientation 
and  scale  of  the  photographs  from  which  the  plate  constants  are  determined. 

A  German  astronomer,   Franz,   established  the  original  eight  funda- 
mental points  in  the  early  1900' s.     Through  the  use  of  five  plates  from  Lick 
Observatory  he  expanded  these  to  a  system  of  150  points.     By  1958,   an 
Austrian  astronomer,   Schrutka-Pechtenstamm,  published  a  revision  of  the 
moon  libration  theory  and  a  recomputation  of  Franz's  150  points.     This 
system  is  considered  the  best  available  and  has  served  as  the  basis  for  later, 
more  densified  systems  [26].  Yet  the  S-P  system  and  others  comparable  to  it 
reflect  the  inaccuracies  inherent  in  the  original  heliometric  observations  as 
well  as  the  additional  inaccuracies  associated  with  the  earth-based  photographic 
process. 

Two  American  government  agencies  have  undertaken  densification  of 


lunar  control.     The  Army"  Map  Service  (now,   Army  Topographic  Command) 
published  AMS-64  consisting  of  256  points.  This  agency  utilized  the  funda- 
mental points  from  the  IAU  Cataloque  of  Blagg  and  Muller  and  plates  from 
the  Lick  Observatory  [8],    In  1966,  AMS  published  the  GROUP  NASA  system 
of  484  points  utilizing  control  points  determined  by  Saunder,    Franz,   and 
Konig  [18].  The  Aeronautical  Chart  and  Information  Center  of  the  U.  S.  Air 
Force  published  another  independent  system  of  196  points  in  1965.     ACIC 
selected  Control  from  the  S-P  system  and  plates  from  the  Pic  du  Midi 
Observatory  in  France  and  the  U.  S.  Navy  Astrometric  Reflector  in  Arizona 
[23].  There  were  large  differences  between  the  systems  of  the  two  agencies  in 
planimetry  (several  kilometers)  and  height.    This  was  emphasized  during  the 
RANGER  probes  to  the  moon  when  elevation  differences  of  approximately  2.5 
kilometers  between  the  AMS/ACIC  systems  and  the  trajectory  computations 
were  noted.     Nevertheless,  the  systems  were  combined  to  form  the  Selenodetic 
Control  System,  DOD-66,  of  734  points  [26][19], 

Two  modern  photographic  methods  are  independent  of  control  establish- 
ed through  heliometric  observations  and  appear  to  be  rather  promising.     The 
Lunar  and  Planetary  Laboratory  at  the  University  of  Arizona  employs  a  pro- 
cedure using  star  trailed  photography  that  was  designed  by  Arthur  [26].    Per- 
haps more  significant  is  a  procedure  contributed  by  Kopal  of  the  University  of 
Manchester.     Moutsoulas  describes  it  as  photographing  a  stellar  field  that 
is  at  the  same  declination  and  hour  angle  that  the  moon  will  attain  at  a  later 
time.    When  the  moon  reaches  the  proper  position,  the  plate  is  reexposed. 
Providing  no  excessive  temperature  changes  take  place  during  the  period  the 
telescope  is  stationary,  the  star  field  provides  the  plate  orientation  and  scale; 
and  the  constants  can  be  used  for  reduction  of  points  on  the  lunar  surface  [24], 
Kopal  states  that  the  achieved  accuracy  is  sufficient  to  determine  the  physical 
librations  of  the  moon  [22].. 

Extensive,  extraterrestrial  photography  was  inaugurated  with  the 
launching  of  the  Lunar  Orbiter  Satellites  during  the  period  August  1^66  and 


August  1167.    The  mission  of  the  first  three  Orbiters  was  primarily 
designed  for  the  selection  of  primary  and  secondary  landing  sites  for  sub- 
sequent  Apollo  missions.     Orbiter  IV  and  V  were  tasked  to  perform  a  broad, 
systematic  survey  of  scientifically  interesting  features  on  the  lunar  surface. 

All  Orbiter  photographic  subsystems  contained  a  medium  resolution 
lens  (focal  length  80  mm)  and  a  high  resolution  lens  (focal  length  610  mm). 
Neither  was  of  photogrammetric  quality.     Calibration  of  the  system,   in 
general,  included  determination  of  the  calibrated  focal  length,   radial  and 
tangential  distortion,  the  principle  point  of  autocollimation  and  the  camera 
format  reference  system  with  respect  to  sawtooth  fiducial s  and  a  preexposed 
reseau  system  on  the  film  (Lunar  Orbiter  I  lacked  these  reseau  marks).     Ad- 
ditional calibration  was  required  to  establish  the  effect  of  an  image  motion 
compensation  system. 

In  operation,  the  film  would  be  clamped  to  the  platen,   and  the  platen 
would  move  in  proportion  to  ground  speed  while  the  shutter  was  open.     The 
film  was  then  processed  by  a  BIMAT  system  which  developed,  fixed,   and 
dried  it.    The  negative  was  then  scanned  by  a  line  scan  tube  in  small  increments 
(2.  67  mm).    This  signal  was  electronically  processed  for  transmission  to 
earth  via  the  spacecrafts'  telemetry  subsystem  as  a  composite  video  signal. 
The  ground  reconstruction  electronics  system  received  the  video  signal  and 
fed  it  to  a  kinescope  tube  from  which  it  was  copied  on  35  mm  television  record- 
ing film.     A  reassembly  printer  utilized  this  record  to  orient  and  project  the 
framelets  on  aerographic  duplicating  film  to  produce  the  finished  product. 

The  photography  collected  from  this  series  eliminated  several  signif- 
icant limitations    attached  to  earth  based  photography;    namely,  the  distortions 
associated  with  atmospheric  refractivity  and  insufficient  scale  for  effective 
resolution.     Further,  it  provided  a  greatly  improved  geometry.     However, 
other  disadvantages  inherent  in  the  total  system  design  requirements  introduced 
distortions  into  the  photography  and  uncertainties  into  the  reduction  procedures. 
Broadly,  the  distortions  were  associated  with  on  board  photographic  processing, 
space  transmission  of  the  video  signal,   and  ultimate  reconstruction  of  the  photo. 
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Reduction  uncertainties  included  the  film  distortion,  but  additionally,  was 
largely  dependent  upon  photo  support  data  which  defined  spacecraft  location 
and  attitude  at  time  of  exposure.     These  were  functions  of  the  orbit  determin- 
ation program  with  its  associated  uncertainties. 

Nevertheless,  despite  the  fundamental  inaccuracies,   ACTC  evaluated 
the  feasibility  of  establishing  a  lunar  geodetic  system  from  Lunar  Orbiter 
photography  and  arrived  at  positive  conclusions  [3].      One  result  was,  A 
Positional  Reference  System  of  Lunar  Features  Determined  From  Lunar 
Orbiter  Photography.     Although  the  original  feasibility  study  encompassed 
only  the  Lunar  Orbiter  IV  Mission  with  its  polar  orbit  and  extensive  coverage, 
it  was  found  that  the  medium  resolution  photography  was  of  particularly  poor 
quality  in  detail  except  near  the  terminator.     The  remainder  was  either 
highly  over  or  under  exposed.     All  photography  possessed  significant  errors 
in  timing,  exposure  orientation,   and  spacecraft  po si tioning [3].       As  a  result, 
photography  from  all  Lunar  Orbiter  missions  was  utilized  in  order  to  achieve 
the  desired  coverage.     However,   Lunar  Orbiter  I  photos  which  lacked  a  pre- 
exposed  reseau  grid  on  the  film  were  employed  only  when  necessary  to  fill  in 
specified  areas.     The  method  used,  broadly,  for  this  control  system  is  best 
described  by  the  author: 

"The  method  consists  of  computing  perspective  projections  [23] 
based  upon  the  orbital  data  for  a  series  of  photographs  that 
are  linked  together  by  common  coverage.    Starting  on  the 
nearside  Lof  the  moonj,  the  projections  were  positioned 
to  agree  with  the  coordinates  of  features  determined 
from  telescope  photography.     [The  ACIC  net  of  196  control 
points,  [23]].    The  link  was  continued  around  the  moon  by 
extending  the  coordinates  of  common  features  from  one 
photograph  to  the 'next.     A  meridional  arc  and  an  equatorial 
arc  were  completed  and  joined  in  the  vicinity  of  the  equator 
and  the  180th  meridian  [27]. 


This  net  produced  (considering  the  extent  of  the  net  and  lack  of  farside 
control)  reasonable  estimated  accuracies  of  1-5.  5-10,   and  10-15  kilo- 
meters, depending  on  the  particular  area  cited  [27].    This  was  achieved 
despite  the  facts  that  control  was  provided  only  on  the  nearside  in  a  coor- 
dinate system  based  on  center  of  figure,   and  the  photography  was  of  variable 
quality  with  all  the  errors  associated  with  its  on  board  processing  and  elec- 
tronic transmission.     Further,  the  exterior  elements  of  camera  orientation 
were  determined  from  spacecraft  telemetry  with  the  associated  orbit  deter- 
mination uncertainties  and  a  coordinate  system  originating  at  the  center  of 
mass. 

A  current  control  net  in  the  process  of  being  established  by  the 
Mapping  Sciences  Branch  of  the  Manned  Spacecraft  Center,   Houston,  Texas, 
is  in  the  imminent  stages  of  completion.     This  net  is  based  on  medium  and 
high  resolution  photography  acquired  solely  by  Lunar  Orbiter  TV.    It  covers 
a  rather  extensive  area  between  ±20°  latitude  and  60°  west  longitude  to  45 
east  longitude  with  the  greatest  concentration  of  control  in  the  Apollo  land- 
ing zone  of  ±5W  latitude  of  the  same  area.     Although  control  points  from 
DOD-66  and  the  ACIC/AMS  nets  are  input  to  the  computational  program,  they 
are  generally  not  used  in  the  adjustment.     They  are  merely  compared  to  the 
control  established  by  Lunar  Orbiter  IV  and  the  root  mean  square  differences 
are  output  in  the  statistical  summary.     Preliminary  results  have  shown  a 
bias  between  the  two  systems  of  approximately  two  kilometers,  but  the  final 
results  have  yet  to  be  published. 

All  of  these  control  systems  are  steps  toward  the  fulfillment  of  the 
essential  requirements  for  the  development  of  geodetic  and  cartographic 
knowledge  of  the  moon  as  outlined  by  the  Falmouth  conference  of  scientists, 
convened  by  the  National  Aeronautics  and  Space  Administration  at  Falmouth, 
Massachusetts  in  July  of  19.65  [12].    Among  these  requirements  are: 

Establish  a  selenodetic  coordinate  system.  .  .  related 
to  the  right  ascension/declination  system. 
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Derive  a  reference  figure  with  respect  to  a  point 
which  is  representative  of  the  moon's  center  of 
mass. 

Establish  a  three-dimensional  geodetic  control 
system.  .  .  in  terms  of  latitude,  longitude,   and 
height  above  the  chosen  reference  figure. 

These  requisites  are  not  only  essential  to  the  expansion  of  geodetic  and 
cartographic  knowledge  of  the  moon,  but  become  fundamental,  base  know- 
ledge for  the  exercise  of  other  disciplines  |~12].    Photogrammetry  has 
demonstrated  uniquely  that  it  provides  the  necessary  capability  to  efficiently 
gather  the  necessary  data  and  to  process  it  into  useful  and  meaningful  infor- 
mation [12]. 

The  following  photogrammetric  procedure  is  submitted  as  a  modest 
contribution  to  the  ever  expanding  numbers  of  methods  designed  to  increase 
man's  knowledge  of  the  lunar  body. 


4.     EXPERIMENTATION 

4. 1    Real  Data  Experiment 

4.1.1  Preliminary 

The  purpose  of  this  demonstration  is  to  describe  in  detail  the 
procedure  utilized  to  transfer  local  or  selenographic  coordinates  from  an 
assumed  or  known  location  to  surrounding  lunar  features  that  are  identi- 
able  in  orbital  photographs.    It  must  be  realised,  however,  that  no  lunar 
surface  photography  has  been  accomplished  with  this  purpose  in  mind.     As 
a  result  several  basic  assumptions  are  employed  and  various  procedures 
inaugurated  that  would  normally  be  unnecessary  were  such  a  mission 
assigned  to  personnel  of  the  APOLLO  series  or  follow-on  series  which  would 
reach  the  lunar  surface. 

4.1.2  Materials 


The  following  materials,  equipment,   and  systems  were  used: 

A.  APOLLO  12  Lunar  Surface  Photographs;    AS12-48-7090,  7091, 
7092;    Magazine  X;    Exposed  by  a  70mm  Hasselblad  camera 
with  focal  plane  reseau  grid.     (Nominal  focal  length,  60mm) 

B.  A.  M.S.,   Lunar  Map,   Surveyor  III  Site;    Scale;    1:2000 
(1st  ed.  ,  Jan  1968) 

C.  Mann  Precision  Comparator,  Type  735  with  Mann  Data  Logger 

D.  IBM  360/75  Computer  System  (OSU  installation) 

The  photographs  identified  in  A.  above  were  the  result  of  an  extensive 
search  of  all  surface  photography  obtained  furing  the  surface  operations  of  Apollo 
Missions  XI  and  XII.    They  were  selected  with  the  following  criteria  in  mind: 

A.  Stereoscopic  coverage 

B.  Maximum  base  between  photographs 

C.  Simultaneous,-  photographic  coverage  of  the  LM,  Surveyor  III; 
and  other  points  on  the  lunar  surface  that  could  be  identified 
from  orbital  photography 
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D.    Exposed  with:  a  calibrated  camera  equipped  with  a  focal  plane, 
reseau  grid. 

These  three  photos  fulfilled  these  requirements  adequately  with  an 
average  base  estimated  to  be  twenty  meters;    the  LM  and  Surveyor  III  were 
imaged  on  each  photo;    three  relatively  well  defined  lunar  features  were 
imaged;    and  a  post  flight  calibration  was  conducted  on  the  two  cameras  em- 
ployed.   Each  camera  was  equipped  with  reseau  grid  at  the  focal  plane. 
Unfortunately,   it  has  not  been  ascertained  which  camera  exposed  these  par- 
ticular plates  [4].     However,  their  calibrated  focal  lengths  of  61.547mm 
(#1016)  and  61.  636mm  (#1002)  determined  at  a  22.  5m  focus  with  black  and 
white  film  (KODAK  S0267)  were  quite  similar  [5]r6].     Neither  camera  had  a 
lens  distortion  pattern  that  would  require  consideration  except  for  the  most 
rigorous  photogrammetric  procedures  [5]r6]. 

For  the  purpose  of  this  demonstration  the  average  focal  length  was 
used  in  calculations.  This  constituted  the  introduction  of  approximately 
±0.07%  error  in  the  focal  length  and  a  proportional  amount  in  the  computations 
associated  with  it.    This  was  considered  insignificant  for  the  purpose  of  the 
real  data  experimentation.     Further  the  reseau  grid  was  assumed  to  be  at 
exactly  spaced  internals  of  10mm,  (4),  and  radial  and  tangential  distortions 
were  neglected  [17]f5][6]. 

4.1.3    Procedure 

Broad  exposure  to  the  many  hundreds  of  photographs  taken  during 
APOLLO  XII  surface  operations  permitted  the  viewer  to  acquire  a  semblance  of 
orientation  in  regard  to  several  features  on  the  lunar  surface.    This  was  not 
facilitated  by  any  documentation  concerning  time,  direction  of  exposure,  orien- 
tation of  the  camera  or  any  other  details  except  in  the  most  general  sense. 
Nevertheless,  this  orientation  permitted  the  selection  of  three  photographs  with 
the  LM,  the  Surveyor  III  and  three  other  photo  identifiable  features  which  could 
be  located  on  the  lunar  maps.     Further,  it  was  confirmed  that  these  features 
could  be  seen  on  available  orbital  photography.     Specifically,  this  was  photo- 
graphy from  Lunar  Orbiters  I  and  III.     APOLLO  XII  orbital  photography  which 
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covered  Surveyor  III  Site  was  taken  at  a  height  of  approximately  60  nautical 
miles  using  a  lens  of  80mm  focal  length.     The  comcomittant  photo  scale 
was  nearly  1:1,400,000.    This  was  entirely  inadequate  for  surface  feature 
identification  within  the  limitations  of  surface  acquired  photography. 

The  lunar  maps  of  Site  Til  were  employed  to  establish  the  coor- 
dinates of  the  five  points  to  be  used.    The  LM  was  plotted  on  Lunar  Map, 
Surveyor  III  Site  (Scale;  1:2000)  from  coordinates  established  on  Lunar 
Surface  Exploration  Map,   LSE  7-6,  Scale  1:5000,  prepared  by  the  U.S. 
Army  Topographic  Command,  1  November  1969.    With  the  top,  center  of  the 
LM  arbitrarily  defining  the  origin  of  a  local  cartesian  coordinate  system  its 
azimuth  from  Surveyor ITI  was  measured  on  map  B  as  301°  30    Q0''0  and 
fixed  to  establish  orientation.     Additionally,  the  distance  between  the  LM  and 
Surveyor  was  measured  and  fixed  at  202.00  meters  to  establish  scale.    The 
local  coordinates  of  the  three  other  points  were  obtained  relative  to  the  LM. 
The  heights  were  determined  relative  to  the  top  center  of  the  LM  by  inter- 
polating between  the  five  meter  supplementary  contour  intervals  provided  on 
the  map.      The  initial  locations  of  all  points  are  summarized  as  follows  (See 
Figures  I,  IA,  and  IT): 

SELENOGRAPHIC   COORDS.      LOCAL  CARTESIAN  COORDS. 
POINT  LATITUDE        LONGITUDE  X  Y  Z 

1  (LM)  3-11-51.6  S  23-23-14.  GW  1000.00  1000.00  100.00 

2  (SURVEYOR)  3-12-04.0  S  23-22-53. 6  W  1172.23  894.46  87.49 

3  (MOUND)  3-11-46.  IS  23-23-20. 3  W  948.00  1045.00  93.96 

4  (LONE  ROCK)  3-11-52.9  S  23-22-58. 8  W  1129.23  988.46  93.96 

5  (CRATER  RK)  3-11-53.5  S  23-22-55. 7  W  1156.23  982.46  91.96 

The  location  of  camera  exposure  stations  provided  a  more  difficult 
problem  since  there  was  no  documentation  in  their  regard.    Therefore,  esti- 
mated positions  had  to  be  determined  from  the  photographs  themselves.    This 
was  accomplished  graphically  be  constructing  a  template  based  on  the  camera 
field  of  view.    With  a  nominal  focal  length  of  60  millimeters  and  usable  camera 
format  of  52  by  52  millimeters  the  angular  field  of  view  was  computed  to  be 
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approximately  46°.     There"  was  an  angular  field  of  9?  2  between  adjacent 
reseau  crosses.    The  template  was  overlayed  on  the  lunar  map  and  adjusted 
until  identifiable  lunar  features  were  in  their  proper  angular  relationship. 
When  the  optimum  fitting  of  the  template  was  achieved,  the  vertex  defined 
approximations  of  the  exposure  station  in  planimetry  (X0,  Y0)  and  the  central 
axis  of  the  template  defined  the  direction  of  the  camera  optical  axis.    This 
provided  an  estimate  for  the  phi  ((p)  rotation.     Exposure  station  height  (Z0) 
was  again  interpolated  from  contour  intervals  modified  by  an  added  1.  37  meters 
based  on  the  assumption  that  the  astronaut  accomplished  the  photography  stand- 
ing with  the  camera  at  mid-chest  level.    Estimates  of  the  omega  (co  )  and 
kappa  (x)  rotations  were  determined  from  the  apparent  depression  angle  of 
the  center  cross  reseau  and  the  comparison  of  a  line  of  horizontal  reseau  marks 
with  the  apparent  lunar  horizon,  respectively.     A  summary  of  the  locations  of 
the  exposure  stations  and  camera  orientation  estimates  are  (See  Figures  I,  IA, 
and  II): 

STATION  SELENOGRAPHTC  COORDS.       LOCAL  CARTESIAN  COORDS. 

(PHOTO#)  LATITUDE         LONGITUDE  Xp  Y0  ZQ 

1  3°12/llf3S      23°22/52?0W      1186.23m       832.46m    94.09m 
(7090) 

2  3  12  09.  OS      23  22  49. 6  W      1206.23  852.96        94.24 
(7091) 

3  3   12  06.7  S      23  22  48.6  W      1214.73  871.46        95.34 

(7092) 

ORIENTATION  (DEGREES/RADIANS) 

X  (£  to 

1  3.50/0.06109*         20.0/0.34907  80.0/1.39626* 

2  3.50/0.06109  42.0/0.73304  80.0/1.39626 

3  3.50/0.06109  60.0/1.04720  80.0/1.39626 
*  A  selected  average  for  the  three  photographs  was  imployed  for  the  xand 

a)  rotations. 

It  became  apparent  during  the  template  fitting  procedure  that  there 
existed  a  definite  possibility  of  a  significant  discrepancy  between  the  location 
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SURVEYOR  III    SITE 


SCALE    1:2000 


MERCATOR   PROJECTION 

STANDARD   PARALLELS   AT    2°30'N   AND    2°30'S   LATITUDES 

CONTOUR    INTERVAL- 10  METERS 

SUPPLEMENTARY   CONTOURS   AT   5    METER    INTERVALS 

CONTOURS   AND   SPOT   ELEVATIONS   ARE    EXPRESSED   AS    RADIUS   VECTORS    IN   METERS   WITH   THE   FIRST   THREE 
DIGITS   OMITTED.    FOR   EXAMPLE:    A   RADIUS   VECTOR   OF    1738250   METERS   IS    DESIGNATED   8250   METERS. 

THE   VERTICAL   AND   HORIZONTAL   CONTROL   NETWORK   ON   THIS   MAP   WAS   ESTABLISHED   BY 
PHOTOGRAMMETRIC   TRIANGULATION    USING   THE    LUNAR   ORBITER  SITE   IP- 7   CONTROL 

RELATIVE    ERRORS    EXPRESSED    IN    METERS   (90    PERCENT   PROBABILITY) 

HORIZONTAL 3  METERS 

VERTICAL     6  METERS 


50 o      meters 


Figure  I 


14 


s 
a 

•l-H 
P>4 


15 


Figure  H.      (Photo  AS12-48-7092) 


16 


plotted  for  the  LM  and  the  position  indicated  by  its  angular  relationship  with 
other  features.    It  appeared  that  the  actual  position  of  the  LM  should  be  some 
25  meters  to  the    NE  of  its  current  position.     However,   since  no  better 
information  on  its  selenographic  coordinates  was  available.     It  was  consid- 
ered to  be  fixed  with  the  qualification  that  this  discrepancy  would  be  investi- 
gated by  varying  the  application  of  constraints  on  it  and  other  points  during 
the  adjustment. 

The  original  intention  was  to  measure  photo  coordinates  on  the 
Zeiss,   Precision  Stereocomparator,   PSK,  with  ancillary  IBM  026  card  punch 
to  facilitate  use  of  the  computer  program  COMCORDCON.    This  program 
converts  comparator  coordinates  to  photo  coordinates  by  an  affine  transfor- 
mation, simultaneously  correcting  for  lens  distortion  and  film  shrinkage 
(See  Appendix  I).    Because  of  the  malfunction  of  this  equipment  the  Mann 
Precision  Comparator  was  utilized.     Unfortunately,  to  simplify  the  obser- 
vation procedure,  each  plate  was  rotated  approximately  30°  to  prevent 
alignment  of  the  measuring  cross  with  the  photographic  reseau  crosses.    This 
prohibited  COMCORDCON  from  properly  identifying  the  four  reseau  marks 
associated  with  each  point  measured  and  correlating  them  to  the  reseau,  photo- 
coordinate  system.    A  simple,  two-point  transformation  routine  was  employed 
to  rotate  the  comparator  coordinate  system  near  enough  to  the  reseau  photo- 
coordinate  system  to  make  the  data  compatible  to  COMCORDCON.    The  out- 
put from  COMCORDCON  was  then  ready  for  input  to  the  BLOCK  TRIANGULA- 
TION  computer  program  (See  Appendix  IT). 

The  following  mean  standard  errors  were  estimated  for  conversion 
to  the  variance-covariance  matrices  for  subsequent  use  in  the  BLOCK 
TRIANGULATION   program  for  weighting: 

Photo  coordinates;         Ox  =  o^=  0. 01  mm 
Exterior  orientation;      aXo  =  by    =  CT.Z  0  =    20.0  m 

CT(0   =  CTu  =    0. 174533  rad  (10°) 
ax   =  0.  08727    rad  (5°) 

Survey  coordinates;        ctx  =  6>  =  az  =  °°to  0.  01  m  (various) 
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4.1.4    Real  Data  Results  " 

In  addition  to  the  variance-eovariance  matrices  postulated  from  the 
standard  errors  of  the  previous  section,  constraints  on  the  survey  coordinates 
of  Points  1  and  2  and  the  elevation  coordinate  of  Point  3  were  imposed  assum- 
ing a  standard  error  of  0.  01  meter.     The  results  of  this  first  adjustment 
were  exceedingly  poor.     Subsequent  adjustments  consisted  of  input  imposing 
constraints  on  combinations  of  Points  1  and  2  and  variable  constraints  and 
relaxations  on  Points  3,4,   and  5.    These  triangulations  either  provided  only 
slightly  improved  results  or  the  adjustment  failed  to  converge  at  all. 

Two  tendencies  were  manifest,  particularly.     Point  1,  the  LM, 
continually  drifted  to  the  lunar  northeast  or  east,  and  there  was  a  constant 
warping  of  the  model  most  evident  in  the  residuals  on  surface  point  elevations, 
the  x  rotation  which  was  constrained  to  5  ,  and  in  u)  which  was  constrained  to 
10°.    When  the  constraints  on  Point  1  were  relaxed,  the  LM  freely  moved 
approximately  48  meters  almost  due  east  of  its  initially  plotted  position. 
The  warping  appeared  to  subside  to  some  extent,  but  further  variations  of  the 
weight  matrices  were  required  to  reduce  the  residuals  on  survey  elevations 
and  the  rotations  associated  with  the  elements  of  exterior  orientation  to  any 
degree  of  realism  or  consistency. 

These  difficulties  were  attributed  to  the  possibly  erroneous  position- 
ing of  Point  1,  the  possible  mis  identification  of  Point  3,   and  the  uncertainties 
associated  with  the  coordinates  of  all  points  that  were  fixed  and  employed  as 
control  for  the  model.    Elevation  differences  were  particularly  noted  to  be 
a  potential  source  of  error  since  the  elevation  differences  among  all  points 
were  relatively  small  and  generally  within  the  predicted  error  of  the  lunar 
map  (6  meters  with  90%  probability).     A  further  complicating  factor  involved 
with  the  uncertainties  in  elevation  determination  and  the  minimal  differences 
was  the  near  coplanearity  of  the  control.    As  explained  by  Smith  [29]  this 
would  manifest  itself  in  the  triangulation  program  as  an  indeterminacy  of  the 
normal  coeffecient  matrix.    Of  possibly  worse  consequence  is  Thompson's [31  ] 
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expansion  of  Smith's  explanation  which  would  indicate,  if  not  indeterminacy, 
then  an  instability  of  the  solution. 

Triangulations  numbered  27  and  29  provided  the  most  consistent 
adjustments  that  could  be  extracted  from  the  real  data.     However,   No.  27 
utilized  Point  3  as  a  control  point  and,  as  a  result,  the  values  involved  must 
be  suspect.    Triangulation  No.   29  utilized  Points  2,  4,   and  5.     Since  these 
points  appeared  to  be  well  identified,  were  in  proximity  to  the  camera 
exposure  stations,   and  had  the  best  known  positional  relationship,  this  trian- 
gulation is  accepted  as  the  most  valid.     Unfortunately,   acceptance  of  trian- 
gulation No.  29  positions  the  LM  at  LAT.   3-11-51.4  S    LONG.  23-23-08. 1W. 
This  reduces  the  distance  between  the  LM  and  Surveyor  III  to  approximately 
163  meters,   and  redefines  the  bearing  to  the  LM  to  about  311°  30  .    This 
possible  redefinition  of  the  scale  and  orientation  of  the  system  effectively 
distorts  the  information  it  produces.     Nevertheless,  the  results  do  have 
value  insofar  as  the  adjustment  retains  consistency  and  merely  lacks  a  valid 
scale  and  orientation.     A  complete  summary  of  the  results  and  the  statistics 
of  these  adjustments  are  provided  on  pages  25  through  36. 

It  can  be  concluded  that  the  possible  gross  uncertainties  of  this 
particular  set  of  real  data  negate  any  reasonable  expectation  of  significant 
results.     However,  the  feasibility  of  employing  real  data  with  proper  control 
seems  to  be  reasonably  apparent. 

4.2  Idealized  Data  Experiment 
4.2.1     Procedure 

In  order  to  provide  a  standard  by  which  one  might  logically  anticipate 
the  predicted  accuracies  of  a  triangulation  program  utilizing  lunar  surface 
photography  with  realistic  control,  an  idealized  model  was  constructed  (Fig- 
ure in).    This  model  presupposes  a  reasonable  capability  of  determining  the 
relative  elevations  of  Point  3  and  the  camera  exposure  stations  with  respect  to 
Points  1  and  2;    and  an  ability  to  make  an  estimate  of  the  xand  od  rotations  of  the 
elements  of  exterior  orientation. 
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I.  (110,110,0) 

1.  (113,132,0) 

II.  (116,110,0) 

2.  (113,132,6) 

3.  (118,132,2) 

4.  (102,160,8) 

" 

5.  (160,300,12) 

Figure  III 
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Additional  conditions  and  parameters  are: 

1.  Points  1  and  2  aligned  with  local  vertical  at  a  fixed  distance 

(In  this  case  6  meters) 

2.  Point  3  at  a  known  elevation  relative  to  Points  1  and  2  and 

at  a  known  distance.     (In  this  case  5  meters) 

3.  A  Hasselblad  camera  (described  previously)  utilized  for  at 

least  two  exposures  providing  a  stereoscopic  pair  at  a 
distance  near  its  22.4  meter  focus. 
The  coordinates  of  Point  1  (assumed  to  be  an  LM  or  other  similar 
landing  vehicle)  are  considered  fixed  and  to  define  the  origin  of  a  local  cartesian 
coordinate  system.     Points  1  and  2  define  the  Z  survey  axis;    and  Point  3  is 
then  defined  to  be  on  a  line  parallel  to  the  X  survey  axis.    The  coordinates  of 
the  control  points  (1,2,  and  3),  the  photoidentifiable  features  (4  and  5),  and  the 
exposure  stations  become: 


POINT 

1.  (Target) 

2.  (Target) 

3.  (Target) 

4.  (Feature) 

5.  (Feature) 
Exposure  Station  1 
Exposure  Station  2 


LOCAL  CARTESIAN  COORDINATES  (meters) 


_x_ 

Y_ 

z_ 

113.00 

132.00 

0.00 

113.00 

132.00 

6.00 

118.00 

132.00 

2.00 

102.00 

160.00 

3.00 

160.  00 

300.00 

12.00 

110.00 

110.00 

0.00 

116.00 

110.00 

0.00 

And  the  rotational  orientations  of  the  cameras  is: 


CAMERA  ORIENTATION  ANGLES  (DEGREES/RADIANS) 

X  (0  to 

Exposure  Station  1     0.  00/0.  0000     0.00/0.0000     90.00/1.5708 

Exposure  Station  2     0.  00/0.  0000     0.  00/0.  0000     90.  00/1.  57  06 
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The  estimated  standard  errors  associated  with  the  various  observations  are: 

Photo  coordinate  :  <jx  =  <7y    =  0.005  millimeters 

Survey  coordinates 

Targets  :  a  x  =  crY    =  ctz    =0.01  meters 

Features:  a  x  =  cty    =  crz    =   °°  meters 

Exposure  Station  :  a  x  =  <7Y    =0°  meters 

az  =   0.01  meters 

ax  =    0.  08727  radians  (5°) 

CT(/D  =  CTC0=:    0.17453    radians  (10°) 

4.2.2    Idealized  Data  "Results 

The  initial  adjustment  of  the  idealized  data  was  a  slight  modification 
from  that  which  is  tabulated.    The  first  triangulation  constrained  the  survey 
coordinates  in  relation  to  the  relative  errors  of  the  lunar  map.    This  was 
assumed  to  provide  standard  errors  of  <j  ^  =  crY     =3  meters  and  crZo   =    6 
meters. 

Although  the  results  of  the  first  adjustment  produced  smaller  standard 
errors  in  the  adjusted  coordinates  of  the  photoidentifiable  features,  the  realism 
of  estimating  the  X0  and  Y0  of  the  camera  exposure  stations  on  the  lunar  sur- 
face to  that  accuracy  appeared  questionable.    On  the  other  hand  the  estimate  of 
elevation  differences  between  the  camera  stations  and  Points  1  and  2  to  a 
reasonable  accuracy  seemed  practicable.    As  a  result,  constraints  on  X<j  and 
Y0  were  removed  and  that  on  Z0  was  strengthened.    These  results  were  predict- 
ably good  and  are  provided  on  pages  37  through  42. 

In  an  effort  to  produce  results  that  might  be  more  indicatory  of 
those  that  could  be  achieved  in  actual  lunar  surface  operations,  the  coordinates 
of  the  surface  features  were  perturbed  within  the  limits  of  the  map  accuracy. 
Expectedly,  the  results  were  identical.    In  a  subsequent  adjustment  the 
constraints  on  the  photo-coordinates  were  relaxed;    that  is,  the  weight  on 
photo-coordinates  was  reduced  from  40,000  to  10,000  (ax,y  =  0.010  millimeters 
vice  0.005  millimeters).     This  caused  a  significant  deviation  of  the  adjusted 
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coordinates  of  the  lunar  features  from  the  known  positions.     In  turn,  the 
constraints  on  the  rotations  of  the  exterior  orientation  elements  were  relaxed, 
and  the  camera  constant  was  perturbed  by  an  additive  0.  050  millimeters. 
The  following  table  provides  these  results  for  comparison.  ■ 

Condition  I:    Coordinates  of  Points  1,2,   and  3  constrained  to  0.01  meter; 
Z0  to  0.  01  meter;    xto  5°;  <p  and  to  to  10  ;  photo-coordinates 
to  5  microns;    and  f  =  60.0  millimeters 
II:     All  of  the  above  except  photo-coordinates  constrained  to  10  microns 
III:     Same  as  II  except  constraints  on  x,  tp  ,   and  co  removed. 
IV:     Same  as  III  except  f  perturbed  (f  =  60.  050  millimeters) 


CONDITION 

ADJUSTED  COORDINATES 

FEATURE  POINT  NUMBER 

X 

y 

z 

KNOWN 

102.00 

160.00 

8.00 

4 
4 

I 

102.009 

159.991 

7.998 

n 

102.010 

159.977 

7.996 

4 

in 

102.010 

159.977 

7.996 

4 

IV 

102.010 

160.001 

7.996 

4 

KNOWN 

160.00 

300.00 

12.00 

5 
5 

i 

159.990 

299.803 

11.986 

n 

159.907 

299.514 

11.967 

5 

m 

159.908 

299.518 

11.968 

5 

IV 

159.908 

299.657 

11.968 

5 

It  can  be  seen  that  the  most  significant  deviation  of  the  adjusted 

coordinates  from  the  known  coordinates  of  the  feature  occurs  as  a  result  of 

relaxing  the  constraints  on  the  photo-coordinates.    This  is  not  unexpected 

since  there  is  a  large  weight  change  involving  the  elements  which  provide 

the  basic  control  for  the  model.     The  only  other  significant  deviation  is 

noted  when  the  focal  length  of  the  camera  is  perturbed  and  this  is  apparently 

confined  to  the  y   survey  coordinate  which  coincides  with  the  rotated  camera 

z  axis. 
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Although  the  scope  of  this  investigation  inhibits  specific  predictions 
of  accuracy,   it  appears  that  with  proper  control  on  Points  1,  2,   and  3  and  the 
Zq  of  the  exposure  stations  a  calibrated  camera  is  capable  of  producing 
positional  accuracies  of  lunar  features  to  several  tenths  of  meters  at  distances 
of  approximately  three-hundred  meters  from  the  control.     The  limited  number 
of  points  negates  any  empiric  estimate  concerning  the  relationship  between 
positional  error  and  distance  from  the  established  control. 
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SURVEYOR  III  SITE 

JOB  NUMBER 

DA  r "E  1*3 "AUG 7~ ; L  97 0 
TIME  15:57:58.2 
""NUMBER  OF  PHOTOS  =~ 
DEGREES  OF  FREEDOM  = 


ADJUSTMENT 
27 


UNIT  STANDARD  ERROR  = 


18 
0.68304D    00 
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? i.O  10    70. 


STO.ERkOR 


RESULTS 
EXTERIOR    ORIENTATION 


"WEIGHTS 


0.1*950-01 

0.00V ~" 


YU    IMEIKKSI 

550. 988 

0.6817U-01 


0.003 


10    (KcTcRS) 

92.701 

0. 1*960-01 

0.003 


KAPPA     (RAO.) 
0.1230960    00 

0.26700-03 
=O7"62OlD=0i  — 
131.312 


PHI      (RAO.) 

0.3669200    00 

0.1 3270-03 


-0.17050-01 
32.828 


OMEGA    (RAO.) 

0.1*26610    01 

0.1**30-03 


-0.30350-01 
"32.828 


?ARIA"iCE/COVAK"lAuC£     MATRIX 
~~ 0.22 3*50-0 3         -G.'87lilD-03           0.1o2o5D-0  3  "' "-0.7*831  0-06  "  ~-0.  5920 10-06" 
~-0. 871510-03            O.*6*670-02         -0.800070-03            0.1625*0-05  0.626960-05 

~~ TJT16T650^O1 '-^HOOtnO'-'OTi 0T22  3CTtT-£n =07T329TJD=7JS ^TfiT  lT9D-0"6"~ 
~-0".7*6E  ID-C6"       ~"0.1h2  5*0-C5         -0. 132900-05 ~~  0.716290-0  7         -0.8  9  3690-08" 
-0.59201Ut06            0.6269t>0-05  "  "-0-851390-06     "    -0.693690-08 —       0.176080-07 
~:Tm"177fcL)-C5 0.529220-05 HT7To9  *  10-05" 


"^0. 117/60-05 
0.529220-05 

"=0T1[S9"*  10^0"5"~ 
"~0. 152130-07 
0.38*120-08 


0.152TTO-57- 


~0"7"3o*12O-0» 


0.2082bO-07 


photo'i^o"; 


~X0~  (METERS)" 
"1192.696 — 


VO"  (METERS)  " 
8o7.*27 


"STO. ERROR 
RESIOUALS" 
VEIGhTS"" — 


0.200^0-01 

■  "ornyjo  02~ 

-         07003 


0.25a"CD="GT' 
-0.1**70   02 
"0.003 


"ZO    (METERS) 
92.515 

—07  5338  0^07" 


KAPPA     (RAO.)  PHI      (RAO.) 

"0.2116130    00  0.67059*0    00 

~~ OT2TTCD-"0"3 


OMEGA    (RAO.) 
0.1195010    01 


"0.L7200   01" 
0.003 


-0.15050    00 
•    "131.312 


"7j776G"30-0* 
0. 62*50-01 
32.828 


0.105*0-03 

0.12*80-02 

32.828 


"0".*"03600-0  3         -07*90360-03" 


-0.*9036D-0"T 
— 07120*60-01" 
-0.37391U-06 


0.6655ltF01T 
-0715*550-03' 

0.16190U-06 


fARIAMCE/COVAKlANCt  MATRIX ~~ 
~  0.120*60-03  -0.373910-06" 
"07T6 1900-06" 


0.635710-06         -0.1001*0-05 


ro7n*"55o-or 

"0.65*100-0*" 
-0.*61220-07~ 


-0.*6l220-07 

0.**5160-07 


-0.356830-06 
0. 132530-06 
-0.658330-08 


0.122690-05 
-0.666280-06 
-0.3*31*0-08 


XI)    (,VET£RS) 
— r2157e*~6" 


STO. ERROR 
"RESIDUALS- 


"'■JETSHTST 


C. 27260    00 

-OVllloO    01 

O.'OOJ" 


YO    (MHTEXS) 
<T6*7tt9  7~ 
0.17*10   00 
_    G.6S630    01 
07003" 


20    (KcTtRSJ 

98.03* 
0.7905U-01 
-0.26990    01 
0;003" 


KAPPA     (RAO.) 
07l9/)6*0~00" 

0.1*150-02 
-0.13610    00 

niT3t2— 


PHI  (RAO.) 
"07965*370    00" 

0.28270-03 
"0.81760-01 
32.828 


OMEGA    (RAO.) 

"O.13*22900r 

0.51970-01 

0.53970-01 


32.828 


VARIA.\CE/COVAKIANCE    matrix 

~D7"7*~2"9"G  5=  01 -C7  *7  G  2  *D~-0"I  0T2 1313  0-0 1        -  07  2 1 02  2  0-~0*" 


"072  9  7  P) 0^0"*  =071X1  63  S=OJ 

~-0.*702*D-0I 07303020-01  ~     -0.136500-01         -0.15*580-0*         -0. 1 305*0-0*   ~  0.8*6980-04 

0.213130-01 -0.136500-01            0.62*8*0-02            0.*690*0-05    "       0.660*00-05  -0.395250-0* 

-U.21022U-U* -0TT5"*53D-T7* 0"."*690*D-05 072C02 00-05         =073*8090-06  -07  l0a**0-06~ 

0.29719D-0*         -0.1305*0-0*           0.660*OD-05         -0.3*8090-06           0.799130-07  -0.3l*>70-07 


-0.131630-03 


0.8*8980-0*        -0.39525D-0* 


-0.105**0-06        -0.31*570-07 


0.270070-06 
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RESULTS 
SURVEY  COORDINATES 


POINT  NO. 

STO._ERROR_ 

RESIOUALS 

WEIGHT 


1019*420 
0.194 OD  00 
0.1942D  02 
0.0 


1029.260 
_0_^1867D  00 
-0.29260  02 
0.0 


100.595 
0.16260-01 


-0.5947D  00 
0.0 


VARIANCE/COVARIANCE  MATRIX 


0. 376200-01    -0.359910-01    -0.129270-02 
-0.35991D-01     0. 348400-01     0.125260-02 


-0.129270-02 


0.125260-02 


0.264440-03 


POINT  NO. 


1172.230 


894.460 


Z 

87.490 


STD._ERR0R 

RESIOUALS 

WEIGHT 


0.38680-02 
0.2  6240-04 
10000.000 


0.49490-02 
(K  72590-05 
10000.000 


0.26280-02 
-0.22930-04 

10000.000 


VARIANCE/COVARIANCE  MATRIX   _^ 

0.149610-04    -0.122270-04     0.233580-05 


0.12227D-04 
0.23358D-05 


0.244880-04    -0.335000-05 
-0.335000-05     0.690860-05 
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POINT  NO.    3         X  Y  2 

948.000         10457000  93.960 


STD.  ERROR        0.6741D-02       0.6647D-02       0. 65610-02 


RESIDUALS         0.5739D-05       0.6810D-05      -0.3354D-04 
WEIGHT  10000.000        10000.000        10000.000 


VARIANCE/COVARIANCE  MATRIX 
"67454370-04    -0.1 731 3D- 05     0. 518780-07 
-0. 17313D-05     0.44188D-04     0.53027D-07 


0. 518780-07     0.53027D-07     0.430440-04 


POINT 

NO. 

4 

X 
1125.3' 

Y 

Z 

tl 

996.982 

90.448 

STD. 

ERROR 

i 

0.99410- 

-01 

0. 

2125D    00 

0. 

,12080- 

-01 

RESIDUALS 

0.3889D 

01 

-0. 

8522D    01 

0, 

35150 

01 

WEIGHT 

0.0 

0.0 

0.0 

\ 

/ARIANCEy 
-02"      -0, 

XOVARIANCE 
,  2 09  5  2 0-6 1~ 
,45159D-0r 

MATRIX 

0.342450- 
-0.74009D- 

-03 
-03 

__. 

0. 

98828D- 

-0. 

20952D- 

-01            0, 

0. 

342450- 

-03         -0, 

,740090- 

-03 

0.145880- 

-03 

29 


POINT  NO. 

5 

X 

Y 

Z 

1156.230 

982. 460 

91.963 

STD.  ERROR 

0. 55500-02 

0.67380-02 

0.55130-02 

RESIDUALS 

-0.35720-04 

-0.74610-05 

0.56370-04 

WEIGHT 

10000.000 

10000.000 

10000.000 

VARIANCE/COVARIANCE    MATRIX 
0V30800D-04      "-0.43  301D-05"    "  -0.382 33 D-06~ 


-0.43301D-05  0.454010-04         -0.200130-06 

-0.382330-06         -0.200130-06         "0.303980-04 
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SURVEYOR  III  SITE 
JOB  NUMBER 

ADJUSTMENT 
29 

DATE  L3  AUG, 
TIME  19:37: 

.  1970 
:25.6 

NUMBER  OF  PHOTOS  = 
DEGREES  QF  FREEDOM  = 

3 

18 

UNIT  STANDARD  ERROR  =   0.68902D  00 


31 


RESULTS 

cxteriur  ume.NrAiio.4 


PHUIU    140.         1 


XU     l.-.C  TLkS) 

lUiTaco 


STO.fcKKUR 


0.  13W0-01 


YO    (MeTERS) 
_o;i6720-or 


20    METERS) 
"91.759" 
—  6".  127  3U-CI 


KAPr-A     (RAO.) 

"~o;i5  3*i80-oi 


PHI      (RAD.) 
~0.  388*060    00 
~  6.13170-03 


OMEGA    (RAO.) 
0.1*57180    01 
"  "0.1**90-03 


Kt'Sl  JUALS 


WEIGHTS 


0.2*30U    01 

oTodF 


-0.2134J    02 
"0.00  3" 


0.2326U    01 
0.003 


0.76*30-01 
~l3i:3l2" 


-0.393*3-01 
"""32.828 


-0.60920-01 
32.828 


VARIANCE/COVAK[Af.CE    MATRIX 
"b".i9S06D-di         -6.063610-03  0.123690-03       "-0.5*0200-06         -o".  364  720-06~      -0.116650-05 

-0.6636UPo"3  oVliTb/Q^OZ         -0. Ts~36*  8  0^0  3  "6TT2  *  SIO-^  7  6T*78~6  30-05        _6~.*680  30-05  ~ 


0. 123690-03 
r0.5*0200-'06~ 
-O.3O4/2U-06 


-0.43e>*S0-03 
If.  12*510-07 
~0.*7863O-05 


0.162110-03 
"-0.799180-06- 
-0.567*50-06 


-0.799180-06 
""0.699500-07* 
-0.125890-0/ 


-0.567*40-06 
~  0.i2569U-0~7~ 
"  0.173420-07 


-0. 165690-05 
0. 101920-07 
0.3*6*20-08 


-0. 116643-05 


0.*6803D-05    -0.165690-05 


0.101920-07 


0.3*6*20-08 


0.2098*0-07 


PHOTO  NO. 


XU  (MEIERS)    YO  (METERS) 


20  (METtRS) 


KAPPA  (RAO.) 


PMI  (RAO.) 


O.*682800-01 


0.7018820  00 


OMEGA  (RAO.) 
IK  1*7  3520  01 


STO. ERROR 


0.200/0-01 


0.2*160-01 


0./2030-02 


0.21820-03 


0.83230-0* 


RESIDUALS 
"htlGMlS""" 


0.11530  02 
OTOOJ" 


-0.13*60    02 
OTOOJ 


0.35960    01 
<H663" 


0.1*260-01 

1317312 


0.31160-01 
32.828 


0.11*90-03 

■67l 7260-01" 

32.828 


0.*02660-03 


~VAR  I  aSiCE/CUVARI  AfcCT  MATRIX 


-0.*5056O-O3 


0.838320-0*         -0.152450-06 


0.693390-06 


-O.*50560-03 


0.583620-03 


0.838320-0*         -0.1001*0-03 


-0.152550-06 


-0. 1S//6U-06 


-0.1001*0-03        -0.197760-06         -0.2779*0-06 


0.516760-0* 


0.1*8770-06 


0.1*8/70-06 


O.*762*0-07 


0.1006*0-06 


-0.103990-05 


0.1216*0-05 


-0.679110-06 


-0.6698*0-08         -0.709520-08 


PHUIU   NO. 


XO    (METERS)  YO    (MEIERS) 


20    (METERS) 


KAPPA     (RAO.) 


Phi     (RAO.) 


OMEGA    (RAO.) 


122*. 119 


862.715 


9*. 512 


-0.3617210-02 


STD. ERROR 


0.272*U    00 


0.14530    00 


0. 57870-01 


0.15900-02 


RESIOuALS 


-0.93B9U    01 


0.6/*5U    01 


0.82270    00 


0.6*910-01 


0.1016510  01 
~ 6.32520-03  ~ 
""0T3069D-0r 


0.1*81900    01 
~  0.7O56D-03 
-0.856*0-01" 


WEIGHTS 


0.003 


131.312 


32.026 


32.828 


VARIANLE/COVARI ANCE    MATRIX 


0.7*1760-01 


-O.*l6U2U-01 


0.152190-01 


0.731960-0* 


O.*09*60-0* 


-O.*16O20-0l 


0.2*12^0-01 


-0.879790-02         -0.8003*0-0*         -0.1512/0-0* 


-0.173610-03 
~0"7lb2  770-03" 


0.152190-01         -O.d/9/90-02 


0.31*920-02 


0.319650-0* 


0.558910-05         -0.39/2*0-0* 


0.731960-0*         -0.8003*0-0* 


0.*09*60-0* 
:07l736iO-:03" 


-0.1512/0-0"* 
"67l02  7  70-"03~ 


0.3l9o40-0*  0.242940-05         -0.37*760-06         -0.4*3120-06 

"o.Via^io-^bs.       -6~.Yi*ib6-d<>        6^105730-06       -6.372210-0/ 

^6"."3972"*D^b*         -0.5*7120-06         -0.372210-07  6".*97900-06" 
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RESULTS 
SURVEY  COORDINATES 


POINT  NO.    I         X     Y  Z 

1048.824         1001.120  112.970 


STD.  ERROR T). J.  373D  00 0 .  1J305D  jOO 0.  2267D-01 

RESIDUALS       _-0.4882D  02      -0.1120D  01      -0.12910  02 
WEIGHT 2lP_  °'CL__ °»°_ 

VARIANCE/COVARIANCE  MATRIX 


0.18839D-01    -0.17755D-01    -0.25725D-02 
-0.17755D-01     0.11028D-01     0. 245220-02 


-0.25725D-02     0.24522D-02     0.51397D-03 


POINT  NO.    2 

X 

Y 

z 

1172.230 

894.460 

87.490 

STO.  ERROR 

0.3915D-02 

0.5040D-02 

0.2643D-02 

RESIDUALS 

0.1376D-05 

-0.1541D-04 

-0.2447D-05 

WEIGHT 

10000.000 

10000.000 

10000.000 

VARIANCE/COVARIANCE  MATRIX 


0 1 . 1 5 3260-04 -0 .12281D-04     0.1501 80-05 

•0.122810-04     0.25407D-04    -0.22320D-05 


0.150180-05    -0.22320D-05     0.69853D-05 
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POINT  NO.    3         X  Y  Z 

"1017.025  993.015  111.559 


STD.  ERROR        0.4686D  00       0.31760  00       0.50790-01 


RESIDUALS        -0.69030  02       0.51980  02      -0.17600  02 


WEIGHT  0.0  0.0  0.0 


VARIANCE/COVARIANCE  MATRIX 


0.21961D  00    -0.148500  00    -0.223130-01 
-0.14850D  00     OT 100860  00    ^0.151190-01 


-0.223130-01     0.151190-01     0.251980-02 


POINT  NO.    4         X                Y  I 

"1129.230          988.460  93.963 

STD.  ERROR        0.6033D-02       0.67030-02  0.5832D-02 

HKFST0U ACS"      ~ 0 . 8 83 7 D-0 4~  "- 0 .  8 05  70- 05  ~  ~07f 31 2 0- 03 

WEIGHT            10000.000        10000.000  10000.000 


VARIANCE/COVARIANCE  MATRIX 
0.364010-04    -0.52303D-05    -0.214820-06 
-0.523030-05     0.44931D-04     0.18533D-06 


-0.214820-06     0.18533D-06     0.34015D-04 
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POINT  NO.    5        X  Y  Z 

11 56.230  982.460  "91.963 


"STO.  ERROR"       0. 55960-02      "0.6790D-62       0.55680-02 


RESIOOALS        -0.9089D-04       0.2998D-04      -0.13040-03 


WEIGHT  10000.000        10000.000        10000.000 


VARIANCE/COVARIANCE  MATRIX 
0.31319D-04    -0.454880-05    -0.322840-06 


-0.454880-05     0.461010-04     0.183230-07 
-0.322840-06     0.183230-07     0.310010-04 
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JOB  NUMBER 


0 


DATE  12  AUG.  1970 
TIM^E  11:39:  3.3 

NUMBER  UF 
DEGREES  OF 


PHOTOS  = 
FREEDOM  = 


10 


UNIT  STANDARD  ERROR  =   0.65281D-02 
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risuuis 
exterior  orientation 


PHOTO  NO.   1 

STD. ERROR 
RESIOUALS 
WEIGHTS 


XO  (METERS) 
109.995 
0.  17**0-0* 
-0.99*60  00 
0.0 


VO  (METcRS) 

110.002 

0.38710-0* 

0.99B5J  00 
0.0 


20  (MEIERS) 

O.COO 

0.15110-0* 

-0.**lOO-0* 

10000.000 


KAPPA  (RAO.) 

-0.1957830-0* 

0.12320-05 

0.19580-0* 

131.312 


PHI   (RAO.) 

-0.2*16350-03 

0.51*10-06 

0.2*160-03 

32.828 


OMEGA  IRAQ.) 

0.1570800  01 

0.52180-06 

0.3*670-05 

32.828 


VARIAHCE/CUVARIANCE  MATRIX 

0.30*080-09  0.370520-09     0.367630-10     0.238570-11 

0.370520-09  0.1*9820-03     0.1*0370-09     0.1837*0-11 

0.3676  30-10  0.1*0  370-09     0.2283*0-09    -0.3*8510-11 

0.238570-11  0.1837*0-11  -0.3*8510-11     0.151880-11 

0.707660-11  0.31*390-11     0.7027*0-12    -0.101370-12 

"0;697200-13~  "  0.1*6610-11  -0.615890-  1 l"~~"  -0.100170-12" 


0.707660-11 
0.31*390-11 
0.7027*0-12 

-0.101370-12 
0.26*290-12 

"0.810650-1*" 


0.697200-13 

0.1*6610-11 

-0. 615890-11 

-0.100170-12 

0.810650-1* 

"0.272270-12 


PHOTO  NO.""  2 


"STD.ERKOR" 
RESIDUALS 
WEIGHTS 


XO  (METERS) 
115.996   " 
"0. 1*1 50-0* 
-0.29960  01 
"   0.0 


YO  (METERS) 

110.000 

0.*1720-0*- 

-0.23**0-03 

0.0 


;0  (METERS) 

-0.000 

"0.15520-0* 

0.39360-0* 

10000.000 


KAPPA  (RAO.) 
-0.15*7180-0* 
— 0.12370-0S" 
0.15*70-0*  ~ 
131.312 


PHI   (RAO. I 
-0.196*960-03 
"0". ;  50090-06" 
~  0.19650-03 
32.828 


OMEGA  (RAO.) 

0.1570800  01 

"0.522*0-06~ 

0.^1630-05 

32.828 


VARIANCE/COVARIANCE  MATRIX 

"0.200110-09-0.7*5**0-10-0.176870-11  ~   0.189810-11  ~  0.577880-11  "  -0.225200-12 

"=0.7*5**0-10     0.17*030-08     0.161900-09    -0. 93*620-12     0.215060-11  0.193790-1 1 

-0.176870-11     0.161900-09  ""  0.2*0780-09     0.537830-11  —  -0. 166670-12  -0.672370-1 1~ 

0.189810-11"  -0.93*620-12""""  0.537830-11"  "  0.152970-11   =0.110780-12  -0.1036*0-12 
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RESULTS 
SURVEY  COORDINATES 


POINT  NO.    IX  Y                 Z 

113.000  132.000           -0.000 

STO.  ERROR        0.8392D-05  0.44970-04  0.03920-05 

RESIDUALS         0.8954D-04  -0. 17210-04  0.67940-05 

WEIGHT            10000.000  10000.000  10000.000 

VARIANCE/COVARIANCE  MATRIX 
0.704300-10     0.42722D-12     0.310480-17 
0. 427220-12     0.20226D-08    -0.823620-15 
0.31048D-17    -0.82362D-15     0.70429D-10 


POINT  NO.   _2         X  Y                Z 

U3.000  132.000            6.000 

STO.  ERROR        0.83920-05  0.4421D-04  0.14530-04 

RESIDUALS        -0.14150-04  0.1749D-04  0. 95270-04 


WEIGHT  10000.000        10000.000        10000.000 


VARIANCE/COVARIANCE  MATRIX 

0.70429D-10     0.421370-12  0.11716D-12 

0. 421370-12     0.195470-08  0.524300-09 

0.117160-12     0.52430D-09  0.211060-09 
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POINT  NO. 


STD.  ERROR 

RESIDUALS 

WEIGHT 


118.000 
0.13000-04 
0.75390-04 
10000.000 


Y 

132.000 
0.44360-04 
0.~2877D-06 

10000.000 


2.000 
0.92830-05 
0.97310-04 
10000.000 


VARIANCE/COVARIANCE  MATRIX 

0.168940-09     0.440220-09  0.39364D-10 

0.440220-09     0.196760-08  0.17593D-O9 

0.393640-10     0.175930-09  0.861660-10 


POINT  NO. 

~~STD~7  ERROR 
RESIOUALS~ 
WEIGHT 


X 

102.009 
0.72980-04 
0.1991D  01 
0.0 


Y 
L59.991 

0.32040-03 
•0.19910    01 
0.0 


7.998 
0.54750-04 
0.15020  01 

0.0 


VARIANCE/COVARIANCE  MATRIX 
0.53267D-08    -0.225560-07    -0.360880-08 
-0.225560-07     0.10263D-06     0.164210-07 
-0.360880-08     0.1642io-07     0.299710-08 
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POINT  "NO.    5         X  Y  Z 

159.990  299.803  11.986 

STO.  ERROR       "0.11460-02  0.4619D-02  0.30070-03 

RESIDUALS"        0.30100  01"  0.21970  01~  -0.9858D  00 


WEIGHT  0.0  0.0  0.0 


VARIANCE/COVARIA.MCE  MATRIX 
0.131320-05     0.528210-05     0.333560-06 
0.528210-05     0.213330-04     0.134720-05 
0.333560-06"    0.134720-05"    0.904040-07 


42 


5.      CONCLUSIONS 

It  is  apparent  from  the  results  of  the  real  data  experiment  that,   in 
general,   the  potential  to  improve  selenodetic  control  by  the  use  of  lunar 
surface  photography  exists  to  a  significant  degree.     Although  the  specific 
results  are  considered  inconclusive  because  of  the  lack  of  any  dependable, 
local  control,   the  experiment  has  emphasized  some  of  the  difficulties 
associated  with  surface  data.      Of  particular  note,   is  the  instability  of  the 
solution  due  to  the  relative  coplanearity  of  the  control  utilized.      This  is 
a  realistic  problem  when  one  considers  that  the  APOLLO  landing  sites  to 
date  have  been  selected  in  the  mare  areas  where  relatively  level  lunar 
terrain  has  been  a  criteria.     It  is  anticipated  that  this  criteria  will 
continue  to  be  considered,  but  perhaps,   to  a  lesser  extent  as  the  expe- 
rience in  lunar  landings  is  increased.      This  does,  nevertheless,   stress 
the  need  for  good  vertical  control,   strongly  constrained,   to  minimize  this 
instability.     Additionally,  the  solution  has  manifested  a  certain  sensitivity 
to  the  rotations  of  the  camera's  elements  of  exterior  orientation.      This 
was  particularly  evident  when  all  elements  of  the  exposure  station  were 
constrained  and  Points  2,   4,   and  5  exercised  total  control  of  the  model. 
The  resulting  adjusted  coordinates  were  realistic  only  for  those  points 
and  stations  within  approximately  fifty  meters  of  the  control  points.     The 
exposure  station  for  photo  AS12-48-7090,   the  most  distant  of  the  exposure 
"stations,   was  almost  two  hundred  meters  from  its  estimated  position  with 
more  than  twenty  times  the  estimated    x  rotation.      Point  3  could  not  even 
be  plotted  on  the  chart. 

On  the  other  hand,   the  idealized  data  and  that  with  perturbations  provides 
some  indication  of  the  kind  of  accuracy  that  may  be  achieved  by  a  reason- 
able effort  to  establish  a  Jocal  network  to  control  the  adjustment  of  more 
distant  features.      Further,   one  may  reasonably  imply  that  an  additional 
input  of  data  from  overhead  photography  (properly  scaled,   if  a  camera 

lens  of  different  focal  length  is  employed)  would  provide  a  material 
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improvement  to  this  adjustment.      Yet,   no  specific  predictions  can  be  offered 
because  of  the  paucity  of  points  and  the  lack  of  suitable  overhead  photo- 
graphy and  information  regarding  the  lunar  conditions  (such  as  surface 
refraction,   etc.).     However,   it  is  justifiable  to  assume  that  the  photo- 
grammetric  errors  associated  with  the  adjusted  local  coordinates  of  lunar 
features  from  surface  and  overhead  photography  would  not  contribute 
materially  to  the  total  error  substantially  attached  to  any  astronomic 
observations. 
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6.      RECOMMENDATIONS 
6. 1     Concept 

The  procedure  to  be  described  is  a  direct  application  of  fundamental 
geodetic  and  photogrammetric  techniques  as  described  in  most  textbooks 
on  the  subjects.  The  unique  aspect  is  that  these  techniques  are  applied  to 
lunar  surface  photography  supplemented  by  orbital  photography.    The 
basic  advantage  of  this  proposal  is  to  establish  control  where  selenodetic 
control  ought  to  be  established. . .  on  the  lunar  surface. 

Although  this  control  will  be  limited  in  coverage,  each  subsequent 
landing  will  provide  a  further  expansion  of  the  control  net  with  an  ever 
increasing  number  of  points  which  can  be  identified  from  orbit  and  to 
which  a  set  of  astronomic  coordinates  originating  with  the  LM  may  be 
associated.    It  is  theorized  that  eventually  a  net  of  sufficient  extent  would 
be  available  to  effectively  control  unmanned,  orbital  photo  missions.    The 
following  procedure  is  offered  to  that  end. 

6.2     Presupposition 

The  current  lunar  landing  vehicles  are  capable  of  obtaining  the 
astronomic  position  of  the  landing  site  from  stellar  observations.    It  is 
presupposed  that  this  capability  will  continue  and  perhaps  improve  in  the 
accuracy  of  determination  as  the  APOLLO  series  progresses.    It  is  further 
assumed  that  an  azimuth  can  be  determined  to  relate  any  local  coordinate 
system  to  the  selenographic  system.    One  method  that  suggests  itself  is  to 
image  a  stellar  field  on  the  lunar  surface  photography  related  to  Universal 
Time  through  spacecraft  time.    This  might  be  accomplished  by  the  use  of 
a  half-circular,  neutral  density  filter  for  the  Hasselblad  camera.    The  top, 
or  clear  half,  would  permit  sufficient  exposure  to  image  the  star  field  while 
the  bottom,  or  tinted  half,  would  inhibit  overexposure  of  the  lunar  surface. 
Time  of  exposure  could  be  recorded  on  a  magnetic  taped  voice  circuit. 
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6.3     Equipment 

The  following  equipment  is  additional  to  what  is  carried  on  the  lunar 
module,   and  serves  only  as  an  example  to  accomplish  the  desired 
procedure. 

1.  A  calibrated  tape  of  approximately  6  meters  that  can  be 
hung  without  interference  from  an  available  or  added 
projection  on  the  LM.      This  tape  would  be  targeted  at 
each  end  with  an  additional  target  whose  position  can  be 
varied  and  its  reading  noted.      A  second,    similar  target- 
for  exposure  station  reference  is  optional.      It  is  visualized 
that  they  would  slide  on  the  tape  with  friction  clamps  to 
maintain  their  position  once  established.      The  lower  end 
should  be  weighted  and  might  have  some  dampening  device 
to  reduce  oscillations. 

2.  Two  lightweight,   variable  height,   telescoping  tripods. 

a.  One  targeted  tripod  would  be  equipped  with  a 
small  leveling  telescope  and  two  calibrated, 
horizontal  spirit  levels.      One  glass  parallel  to 
the  telescope  optical  axis,   the  other  normal  to  it. 
A  plumb  bob  or  optical  plumb  is  necessary, 

b.  The  second  tripod  would  provide  an  attachment 
for  the  Hasselblad     camera  with  similarly 
oriented  spirit  levels. 

3.  A  calibrated  tape  of  convenient  length  (perhaps  up  to  20 
meters)  with  staking  rings  at  each  end  and  a  tension 
spring  with  scale  at  one  end. 

"    6.4     Procedure 


At  any  specified  time  during  lunar  surface  excursions,   the  astronauts 

would  carry  out  the  following  procedure: 
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1.  The  vertical  tape  would  be  hung  from  the  LM. 

2.  Within  20  meters  of  the  LM  on  reasonably  level  terrain 
that  would  include  a  background  with  a  maximum  number 
of  discrete  features,   set  up  the  level  telescope  in  such  a 
manner  that  it  is  level  and  its  field  encompasses  a  portion 
of  the  vertical  tape.      Position  one  of  the  adjustable  targets 
so  that  it  is  centered  on  the  telescope  cross  hairs. 

3.  Lay  out  the  20  meter  tape  from  the  base  of  the  vertical 
tape  to  a  position  below  the  plumb  of  the  level  telescope. 
The  tape  may  be  staked  in  position  with  a  predetermined 
amount  of  tension  indicated. 

4.  Position  the  camera  at  its  first  station  such  that  its  optical 
axis  is  perpendicular  to  the  vertical  tape,   though  not 
necessarily  in  the  same  plane.      Include  in  its  field  of  view, 
the  vertical  tape,   the  targeted  level  telescope  and  the  desired, 
discrete  lunar  features. 

5.  When  the  oscillations  of  the  vertical  tape  are  minimal, 
expose  the  plate  and  record: 

a.  The  reading  on  the  20  meter  tape  below  the 
vertical  tape. 

b.  All  spirit  level  bubble  positions. 

c.  The  reading  on  the  20  meter  tape  below  the 
level  telescope  plumb. 

d.  The  readings  of  the  variable  target(s)  on  the 
vertical  tape.      (All  readings  could  be  voice 
recorded  on  tape.) 

For  subsequent  exposures,    it  would  only  be  necessary  to  reposition 
the  camera  to  obtain  a  stereoscopic  pair,    possibly  readjust  the  optional 
variable  target  (if  used)  and  to  record  the  readings  already  mentioned, 
(See  Figure  ni). 


The  resulting  models  would  be  similar  to  the  idealized  model  described. 
The  vertical  tape  would  provide  the  scale  to  an  estimated  accuracy  of  a 
few  millimeters  and  define  the  local  vertical  as  the  local     Z    survey  axis; 
its  lower  target  would  establish  the  origin  of  the  local  coordinate  system. 
The  targeted  level  telescope,   the  position  of  the  variable  target  on  the 
vertical  tape,  and  the  measured  distance  to  an  estimated  accuracy  of  .01 
meter,    could  be  computationally  corrected  to  define  a  line  parallel  to  the 
local    X    survey  axis.      The  local    Y    survey  axis  would  then  be  defined. 
Reduction  of  the  recorded  readings  would  give  the  survey  coordinates  of 
the  level  telescope  target  and  the    Z    survey  coordinate  of  the  camera 
stations  to  an  estimated  accuracy  of  0.01  meters.      The    x  and    oo 
rotations  would  be  estimated  to  be  near  zero  based  on  the  camera 
leveling  results.      The    cp  rotation  would  be  estimated  by  its  relation  to 
the  defined    YZ    survey  plane.     Approximate  positions  of  the  lunar  surface 
features  can  be  scaled  from  a  convenient  lunar  map. 

After  preprocessing  the  necessary  information  and  providing  the  photo 
coordinates  from  COMCORDON  to  the  BLOCK  TRIANGULATION  PROGRAM, 
the  resulting  adjustment  would  photogrammetrically  relate  all  discretely 
imaged  lunar  surface  features  to  the  position  of  the  LM  in  a  local  cartesian 
coordinates  system. 

Extending  this  with  an  azimuth  and  the  astronomic  position  of  the  LM, 
this  adjustment,   with  a  simple  coordinate  system  transformation  program, 
would  provide  selenographic  coordinates  and  relative  elevations  of  lunar 
features  that  could  be  related  to  current  and  future  orbital  photography. 

It  is  acknowledged  that  the  foregoing  method  is  neither  the  most 
simplified  nor  the  most  sophisticated  that  could  be  employed.      However, 
it  does  serve  to  emphasize  the  fundamental  requirements  of  the  system; 
that  is,   the  establishment  of  adequate  scale  and  orientation  and  the 
application  of  sufficient  constraints  to  obviate  coplanearity  of  the  model. 

It  is,   therefore,   suggested  that  the  previous  procedure,   or  one 
fulfilling  the  same  basic  criteria,  be  considered  for  adoption.      It  is 
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firmly  believed  that  its  implementation  would  be  the  beginning  of  an 
improved  selenodetic  control  network. 
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APPENDIX  I 


COMCORDCON 


The  purpose  of  COMCORDCON,   an  acronym  for  Computor  Coor- 
dinate Conversion  Program,  is,   as  its  name  implies,  to  convert  stereoscopic 
comparator  observations  of  object  space  points  on  a  photographic  plate  to 
photographic  coordinates  in  a  photo-coordinate  system.    The  resulting  measure- 
ments,  sequentially  corrected  for  film  /  emulsion  shrinkage  and  radial  lens 
distortion,   are  then  in  the  proper  preprocessed  form  as  partial  input  to  a 
simultaneous  block  triangulation  program  to  be  described. 

To  achieve  this  objective  COMCORDCON  utilizes  in  this  case  negatives 
or  diapositive  plates  from  a  calibrated  HASSELBLAD  camera  equipped  with  a 
focal  plane  reseau  grid  and  nominal  60.  0  millimeter  focal  length  lens;    the 
Zeiss,    PSK,    Precision  Stereocomparator  augmented  by  an  ancillary  IBM 
026  card  punch;  and  the  OSU,    IBM  360/75  computer  using  a  FORTRAN 
IV  G,    Level  18   Compiler.      Comparable  systems  may  be  substituted  with 
necessary  modification  of  data  format  to  provide  compatibility. 

Input  consists  of  the  number  and  interval  spacing  of  the  reseau  marks 
on  the  focal  plane  grid;    the  type  and  identification  of  the  plates;    lens  distor- 
tion values;    the  coordinates  of  the  center  crosses,  the  object  space  points, 
and  four  neighboring  reseau  marks  (preferably  surrounding  the  object  space 
point);    and  other  job  keeping  information.    This  information  is  then  provided 
to  the  SUBROUTINES  TRANCP  and  RADIS  1. 

TRANCP  performs  an  affine  transformation  by  solving  the  trans- 
formation parameters  for  each  set  of  four  observed  reseau  marks  associated 
with  an  image  point  by  a  least  squares  adjustment.    These  transformation 
parameters  are  then  employed  to  obtain  the  photo  coordinates  of  the  image 
point,  corrected  for  film  shrinkage,  in  the  reseau  coordinate  system  of  the 
photograph.    The  mathematical  formulation  for  this  transformation  is,   for 
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example: 

Assumptions  1.    Peseau  interval  is  exactly  10  mm. 

2.     No  correlation  between  observations  of  unique  points  and  no 

correlation  between  measurements  of  x  and  y  coordinates 

in  an  image  point  observation. 


Y  1 


COMPARATOR 


(FOR  NEGATIVES) 


X 


X=A0  +  A1X'  +  A2Y' 
Y  =  BO  +  BIX'  +  B2Y' 
Fx=    X-A0-A1X'-A2Y'  =  0 

Fy=    Y-B0-B1X'-B2Y'  =  0 


BX  = 


B  Y  = 


±1* 

a  a 


5FV 
SB 


-1 

-xl 

-Yl 

-1 

-xh 

-Yi 

-1 

-xj 

-Yl 

*3 

-1 

-xl 

-Yl 

*4 

E  X  = 


■M 


OBS 


E  Y 


[-3 


f- 

F 

Fxs 

F*3 

F*4 

L.           J 

OBS 


W  =  I;      C  =  BTB;        UX  =  BTEX;  UY  =  BTEY 
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PARX  =  -C*UX 
PARY  =  -C'UY 


A  =  A  +  PARX 
B  =  B  +   PARY 


The  image  photo  coordinates  corrected  for  film  shrinkage  in  the  reseau 
coordinate  system  are  then: 


v<        A2  (Y  -  BO)   -  B2  (X  -  AO) 


A2  Bl      -  Al  B2  [Where   X  and   Y  are 

the  average  of  four 
.        B1(X  -  AO)    -    A1(Y  -BO)  .  ..  , 

Y    =   i ' 1 '  observations.] 

A2B1       -  A1B2 


Assumption  No.   1  is  not  necessary  since  the  reseau  marks  can  be  measured 
individually  and  their  true  coordinates  read  into  the  program  [17].  The 
standard  error  of  their  position  can  be  determined  and  propagated  through  the 
complete  system  if  desired. 

Assumption  No.  2,  though  not  entirely  true,  can  be  minimized  by  careful 
attention  to  proper  observation  techniques. 

Both  assumptions  are  made  for  the  sake  of  simplicity. 

SUBROUTINE  RADIS  1  employs  the  current  image  coordinates  and  applies  a 
suitable  correction  based  on  camera  calibration  data  for  radial  distortion  as  a 
function  of  the  radial  distance  from  the  center  cross.    It  assumes  a  symmetric 
radial  distortion  and  that  the  center  cross  is  coincident  with  the  principle 
point  on  the  plate  [17].  Again  those  assumptions  are  made  for  the  sake  of  sim- 
plicity.   It  would  be  entirely  reasonable  to  integrate  into  RADIS  1  a  correction 
for  the  deviation  of  the  center  cross  from  the  principle  point,   a  model  for 
non-symmetric  radial  distortion  or  a  thin  prism  model  for  correction  of 
decentering  and  radial  distortion  [11]   [10]. 
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COMCOPDCONsupplements 

P-oviding  the  unit  standard  *•  -  oorrectea  photo  eoor(Jinates  by 

-—on  md  fte  st£mdard  o ;- ;;— — the  ^  trans. 

object  space  point  observed.  ^  Y  °hsm 


ons  on  each 
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COMCORDCOM 
LIST  OF  PERTINENT  VARIABLES  AND  INPUT  DATA 


Card  No. 

Variable 

Purpose 

CR 

True  Reseau  Coordinates 

A,  B 

Transformation  Parameters 

1 

INFO  (1) 

Job  IdentificationNumber 

INFO  (2) 

Reseau  Interval  (in  microns) 

2 

CC  (1) 

X  coordinate  of  left  plate 
center  cross  (mm) 

CC  (2) 

Y  coordinate  of  left  plate 
center  cross  (mm) 

CC  (3) 

X  coordinate  of  right  plate 
center  cross  (mm) 

CC  (4) 

Y  coordinate  of  right  plate 
center  cross  (mm) 

COND 

Defines  the  plate  (s)  being 

TYPE 

INFO  (4) 

INFO  (5) 

CCO  (1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

(11) 
(12) 


used,   (RIGHT,  LEFT,  or  BOTH) 

Defines  the  type  of  plate, 
(POSITIVE  OR  NEGATIVE) 

Left  plate  photo-identification  number 

Right  plate  photo-identification  number 

Point  Number 

Point  X  Coordinate,   Left  plate 

Point  Y  Coordinate,   Left  plate 

Point  X  Coordinate,  Right  plate 

Point  Y  Coordinate,  Right  plate 

Blank 

Point  Number 

X  Coordinate,  Upper  Left  Reseau(#l)  Left  plate 

Y  Coordinate,  Upper  Left  T*eseau(#l)  Left  plate 
X  Coordinate,  Upper  Left  Reseau(#l)    Right  plate 

Y  Coordinate,  Upper  Left  Reseau(#l)  Right  plate 

Blank 

59 


Card  No. 

4  Repeat  Observations  of  Card  #3  but  Observing  Upper  Fight  Peseau  (#2) 

5  Repeat  Observations  of  Card  #3  but  Observing  Lower  Right  Reseau  (#3) 

6  Repeat  Observations  of  Card  #3  but  Observing  Lower  Left  Peseau  (#4)  - 

The  set  of  cards  (3, 4,  5,  6)  is  repetitive  by  N  number  points. 
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COMCORDCON 

lyPHCll     •<(?/»!.♦>)      IA-H.I-/I 

III    t/l     in       .    n..         _...  "   *  *     » 


GOIyJ    , 

C0o>  ^J!;f;o'^0«''».ffJB,INCF<(4J 

I;?" uii,j 


0006 
CO';/ 


1 >=-  >0.c»Fl.LAf  I J J«lc.C 
CQii —  ,      lb    "'"•''pltOII     I  INUJI  1  )  .1-1     v> 


000  J 
0010 


?0K^z:'°H^iS«:«:---:-^:^::^«-'-H 


CO  14, 

"016  loco 


UQ     JO" 1=1, a 


-»     ;y°°  ss*  irAi^rj^wKMsaj.^^^: 


C.  IFbl     FLK     pnl;    l;F    C4,A 

£. !§sj_J.vi»fc  cf  iNHtr 


COla 


C 


co?o 

C0?1 

_00«^ 

C023 


35  I'Mfgi  n=f,pci.Nr-i  — - 

■a — r^0?!&^ 

g""  "iKtilHj^^^M'J" 

^°/9  <.0    fc(JC  =  l 

JJJJ  ...    _.  U>    TO    J5 

"«  ioc  j,rcj>      - 

™« ENO 
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IRANCP 

0001  SUUhOU  I  INh     TrtANv.PI  INFI."  ,CCL>  »CP  ,UR  ,CC  ,C('Hl)»  TYPt) 

000?  IMPLICIT  RFAl»o  (A-M,(,-2) 

000*     '•'     ""  '     KfALci    Cl).\i;tUI0Hf /•«  JGt-T       "    ,'7Yl»UTH>*3CIH             '•/VLf-FI/'LrFI            '  •/ 

000'.  LA  I  A     I'l.l/  •  I'M1.                  •/ 

00o>  DIMENSION    INFUCj)  ,cCC(t>0,<i,l?),Cf  ('iC,  12)  ,aC,3)  ,fcXC)  ,EY(A)  ,bT(  J,<i 

*)»UXI  n  iLYIU.Cd,  1)  .Cim.WliCU'.  ).lflJI<M'.,4).l,AUX(3),PAKY(3> 
0006  INTfcCfcrt    lUPCS    l*)/4,4,  1,2,1/ 
_C_ 

C  UCICKMNE  "lht     PLAILIS)     ftHllSE     CC.f.KUlNA  I»  S     AUt      IC     Ul:      iKANSfCrtfLU 

(.  IP     LUNO=iUJTH,     Ull  Tn    K1«JHI    AND    LfcF  I    CCL*  j  J  .*  A  !  i- S     *ILL    Uc     TWANSFCKKFO 

C"  -  IF    LL,.JL  =  Lrt-r,    CrtLY    IHF    LEU    I'LAlf    CLOKuINAlES    WILL    be     TKANSFURMEu 

_C  IF    CUNC=KIGHI,    UNLY    THE    KIGHI    PLATE    CCUKDINAICS    WILL    BE 

C  IKANSrUKPtO 

C_      CCNC  =  LEFT_IS_  ASSUMED    AS_    ThE_  DEFAULT     VALUE  . 

C  . 

C007  INIT*0 

0008  IF  tCHNU.CC.KICihf  )     INII=2 

C 

C  JETFRMINE  IF  rhE  PLAIE  IS  POSITIVE  CR  NEGATIVE 

C  TYPE=NEG  IS  ASSUMED  AS  THE  DEFAULT  VALUE 


OCO'J IDTYPI.  =  2 

0010  1FITYPE.EU.P0S>  "ICIYPE=1 

COU _.\PCINT  =  iNFGt  J) 

0012  R'l'="li\F  C '{  2  )  /  to*  •  3 

0013  J?C1NT-=C 


C014  ICUMI'  =  0 

COlS N  =  l 

C016    ""  10    JPUlNI=JPUlNl"*l  '  •'  •"" 

COll 15    IFUPUINf  .UT.N-i'UlNT)    GO     TO    20C 

C 

C_  HSIAHLlSH    A     THKPC«ArtY    AKKAY    CUNT A  IN ING    THC    COUKDINATCS    0F_THE 

C"  "PuIM     10    bE     IKA.MSFCKtfEC    "IN  10     IHfc   "kEScmU    SVSTLM 

_C  

C018  ""   '"'  CC    20    l«l,V 

0019     TFCKMC I, l)=CCU< JPUINI, I,  INIT»2) 

0020  '  IFCKKl  I  ,2)=Cv.L{  JPUlNl  ,  I  ,  INIT*3) 

_0021 TFLKM1  I,3)«CCU(JPQINT,  I ,  INI  I+jb) 

C022  2C  IFUilMl"l,^)=CL,"U(  JPU  INI  ,"f  i  IN  I  f*9  ) 

0023 __ Kl=INir*l 

0024  K2=1MI*2 

C 

0  '"  COMPUTE     THE    KFSFAU   "lUENT  IFICATICU    NOM.KER    AND    CSTARLiSH    THE    CK    MATRIX 

C 


C02i  AC  =  CCUl) 

C026  UC-CCIK2) 

C02/  Al=l.O 

0028  _^               A2=C.O 

0029""  IU  =0.0 

COiO  rt2=l.O 


C031  SUMX=0.0 

QUi?  •       m  SLMT*0.0 

00  i  J  i'J    uC     100     J«l,', 

003A IiH  =  I?C.C»TFi:RKU,3)-CCUm/Rl»l.*  _  __    :.._.... 

003$  U)2=|?0'.U-[Fi]KH(.l,'tI«CC(K2)  1/iUO.O  "" 

C , 

g  if  \tctss"AkV,  "cuirvtKr   to  thf  positive  sySIcm 

c 

0036  IF{  101  YPC.EO.  I)     UU» I0PUS1 101 > 

C03f  •  CCKA^CKt 102,101, I )»(-l.0>»*I0TYPfc _ 

c 

C  C-CMPU I  F  VALUES  OF  X,Y,  ASP  ESTAliL  1 SH  THE  EX  AND  CY  KATHICES 

C  ASSUME  ThQ.    FOLLOWING  INITIAL  CONfXft  IUNS 

C  _  A0  =  XCC  '  

C Of)=YCC  ""   "~ 

C  Al  =  ti2=l.U  _  _  ■ 

"C~  A2=Kl=C.O 

C  ___._  
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~0d  in  x"=  AC+ n  *c*orx*  aT*Cr"  n  02 1 1  ffl  ,  ?T 

C039  _   '    Y  =  QC*Ul*CnRXtU2*lK(  1U2_|1UI»?) 

C040  EX(  J)  =  TFG<MJ,3)-X  " 

00<>1  tYI  JJ  =  TFUKMJ,4)-Y 

oo<,2  '  iFncoKP.L-c.i)  cu  roico" 


C"  ~E*S f  a'dlTST  The "  b  >.  aTk  l X 
c 

CiHi  ii(J,ll=-1.0'" 

C0<.<.  U(  J,2)  -CCKX»(-l.O) 

CO<;5""    '"  b(  J,JI=CH<  102,  IU1,2)*(-1.0") 


c"  TownrTife""^  "tfF  rnriiB jL'<rrTf:A"(re~x"(rcr«o'[NA"Tis 

.  c  

C(H6    "  ~~SU,VX  =  SUMX»TFCXM  J,  U 

0047  SU.VY  =  SUMY*TFUKM(J,2) 

■  C  "  ...  - — 

C  TEST  FCK  FOURTH  RESEAU  POINT.   IF  TRUE  CONTINUE.   OTHERWISE, 


C  KEiukN  fcTft  NcxTTTiETtACj  PUfRT 

c_ 

00<.0     "  100    CUNIINUE' 

00*9 _  IF(  ICUMP.EC.l)    GO    TO    125 

"C  '."" 

C COMMUTE  THE  AVEKAOE  UF  THE  OBJECT  SPACE  IPAGS  COORDINATES 

C 
0050  XAVG=SUMX/4.0_ 

0051"  YAVC»SUKYA.d 

C 

C  TRAfsSPCSF    KAfKl'X    H 

C 


C0s2  CALL    UGMRAfh.bl  ,4.3) 
C 

C  (JUIA14     l»:f     IJX     A.*U    UY     f  A  I  ■<  I L  F  S  ,     fHlKC     UX-»  ll  IK  AnS  »CX    A*J    U Y *U I R ANi«C Y 

C_ 

covi  "CAUDCPPi<nutr,Ex,iix,"i,4,  i) 

_C0j>4 CALLOGMPkD ( III  .L-YfUYt  if  4?  t J 

c 

C  COMPUTE  C  MATRIX                                 ._„_..  _ 

C 

0055  •  CALLUCMPKOtBI.b.C,  J  ,  <. ,  3  ) 
C 

C  COMPUTE  IWERSG  OF  C  MaTRU 


C0>6  CALIUMINV(C,3,C,PARX,PAKY> 
X 

C  COMPUTE  IhE  TRANSFORMATION  PARAMETERS 

C'  "FORK  HAR  MATRIX,  hhERE  PAR=- I ,C»C INVRS^U" 


C 


00)7  00  JO  1=1,3 

0058 ;   00  bO  J=l,3 

C059    '  "      80  CI  I  ,J)  =  (-l.C)»CU, J) 

C060 CALLUCiMP<C(C.UX,PA*X,  3,3,  I) 

C061  "            '  CALLUjMPUUlC,UY,PA«Y,"3.3,  I)" 
0Co2 ; »\C=HARX  (  I  )^A() ^_ 

~0C6J  Al"=PARX(2)  »aT 

OOoA A2=PA«X{  1)  *A2  

0065  ""  BC=PARY  I  l»*BO 

00o6  _     BL=PARY(2)*Bl  _ 

"C067"  U2-PARYI il  »H2 


C  tOTPini  "cOURUlNAlCS  IN  RbStAU  SYSTEM  UPR  1ME  .  YPR  tME  > 

_„__     C  _  

0C68  "'  UENCM=A2*Hl-Al*b2 

C0o9 TERKl  =  YAVG-00  _. _...... 

CO/0  '"  '     ~"~  Tb<<M>  =  XAVU-AU 

C071 XPK  irt  =  !  A2*TL'i<t'l-82»TFRM2>/DFNCM 

"C072  YPR lKt« (  o  I*  I  EJ<K2-A  I  ♦  TERW  I  )7ut  NUM 

C073 CALL    KACIS1 (XPRlMt.YPlUMb.XP.YP.OISTOR) 

00/<.  "     '  *.i  =  b*lSlt/2 
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c 

C                L'STAOLl'Sh    CP(  J,K)">  AlKlX,     WHERE    J  =  J-lM    PulM,     ANO    K*K-TH    TERM 
C  FOK    K=  I     CK     7  PlililO     IDCNI  IF  ICAT  ICN    NL-KHER 


C  FOR  K  =  '2  OR    8  POINT     IDENTIFICATION    NCKbCR 

C  _FOR  K  =  3  OR    9_  X    PKlMF.    V.ITH    THE    RAOIAL    OISTCRTIUN   CCRRECT1CN 

"C*~  ""'"  APPLIED 

C  FCR  K=<i  CK    10  Y    PRIKE    With    THE    RADIAL    OISTCRTION   CCRRECTICN 


C  APt'LlFU 

C  FOK    K  =  'j    CK     11  SIANUARIJ    UN  II    EKRCR    FOK     IHE    FCUK    RESEAOS    tNLY 


C  FOR    K.  =  6    CR     12  SIA.NUARL7    ERRCR    CF    CNfc    ObSERVATION    ON    TmE    CUJECT 

'       (.  SPACE    POINT 

C 

00  75 CP(  JPIJIN1  ,  UK  5)=  INFO!  A*  INI  T/2) 

00/6  UMJiMINri?»Kn«CCCtJ»»ClNTiirn 

CO// CP(  JPUlM^^iOMXP 

00/8  CP(  JPuT'M  ,4»fci')*YP 

00 1'i  ICCM'M 

OOttO  OU    TO    25 

C_ 

"C    "  CCKPuTt  STANOARO  t'RKCft  -  POINT  BY  PCINT 
C 


OubT  123- TEhi'-O.C 

0082  IEKP1=0.0 

COtti  "  00    ISO    l«H 

COS*  UNPUrEePi  +  nFCMfl  l,l)-XAVG)"2t(TF0RM<  t  ,2)-YAVG)**2 

OOo*;  IOC    TL.it'-  TLK«'«L<(  I  )**<»«fcY(  I  )**2 

0066  tP(JPUlNT,5*K3)=LSCXTUEMP/2l 


C0b7  C  P  (JPG  I  'NT  76  "•»'*  D  =0  SC  kTITEmPT  /<.  8 . 0  J 

C088 'lCOri»»C  __  - 

C0tt9    "  ""  lHCUMO.Mf.BCTH)    GC    TO    10  ™ 

0090 JPC"lNT=JPUl.NT»  INIT/2 

0091  " N  =  N»l 

0092 INI  T=IMI*(-t  J«*N»2  __, 

"0093  '   oO  Tu  !"> 

0094 200  hKITl  (6.20CA)  INFOtll 

0095  2C0A  FOKHATI  '  l«,T52,»JOL  NUMhFR', 16.//T17, 'PHGTC  COORDINATES  CCRRECTEO  '" 

»FCR  ILNS  AND  FILM  OlSICKIlONS  ULISS  KKK-AK- I  5/23   NO.  21  197)',// 

♦/fbC'OMT  S(A<*UAi<U  EKRUK  (KM)  • , IbO, 'SIANLARU  ERROR  OF  M.fcAN  OF  X  A 
+.V0    Y'  t/T  10,  •kHOTC  t  T16i  'POINT  •,  T2tt,  'X  <MM)«^T39,«Y  (KM)  '  ,T<,9,  'AFTE 

*k  Ai  F I  nf   I  k a.»s»  fiu'MA  f Ton  ■",  I  bTTTJn   i  iitT'Ca JECT   5PaC£  >0INT~iW*Tn71 

C096      _j IF JCONU.EC.RIOM  )    GO    TC    225 

C097  "  V.RITEI6.2CC5)     (  (  CP  (  1  I  ,  J  J  )  ,  JJ» 1,6)  ,  I  I*  1  .NPCINI  i 

C098 00    210     l=l,NPOIiU  

C099        "  ""'  IPHCIO=CP( 1.1 ) 

0100 IPCINI =CP( I ,2) 

"010 1  ItMP  =  0SC-<T(Cijri,i>»»2*CPTl,6)*»2) 

0102  21C  PUNCH     JCCO,  IPMIfO,  1P0INT,  (CP  (  I, J). JO. 4)  ,TFMP 

0103  "      ~3CCC    FCKl'AT  (2i'>,  3F  10. A) 

CIO* 225    IFILONO.EO.LCFU    GC    10    250 

0105  "  "  wRlTLUwOO'J)     (  1CP(  lI,JJ).JJ  =  7,l2),U=»l,NPCINT) 

0106  uO  Z  i'i     I'ltNPOlNI 


0107  IPI-.ClU  =  CP(  1,7) 

0108  1PCINT=CP(  1,8)  

0109  '"  ~    '  "235  PONCH  3CC0,  IPK.CTO  ,  I  PO  INT  ,  (CP  (  I  ,"J  )'.  J»"9  ,10  >  ,  TEMP 

0110        200S i_  FORMAT ( T 10 .F6.C, T 18 ,F6. 0.T26.  F9.<i ,  1  37  ,F9. A,  153,615.5,  T 67, EI  5,5, /(T 

250RETURN 
END 
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•4AUIS1 


COOl 

CO02 

"OOUi" 

ccc<.~ 

COO'j 

0006' 
ecu/ 
cooa" 
cooo 

COIO 

con 

CG12 
C013 


SUOKUur  I.»C  I<U;IS1(XP,YP,XP2,YP2,DISTGR> 
IMPLICIT  «EAL*H  IA-H.G-Z) 


Ki-ai      ulc.  .w<,(|<,  )/O.C,U. 0,0. O.O. C.G.G. O.O.O.O.O. Q, O.G.O.O, u. 0,0. 0,0. 

10  .0  .0,0. L, 0.0/ 

rtAO=UitftI IXP**2»YP**2V 
_ If (HAO.Nt.C.C J    GG    TU    IC 

oisroH=c.o 

RETURN 


1C    K3KAL>y\10.0f  1.0 

.  UISTUi^UlSAVG(K)/lC*»3_ 

xp2=xp-cistg*»xp/.<au "" 

YP2=Yf-CISI(!R*YP/RA0 

RETGRN 
ENC 
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*  **     SAMPLE    DATA     COKCORDCON     **  * 


1 

left 

\'rG 

7r  a^ 
C  1 

D  CO  C"  V. 
-  23  7^7 

00 

1 

CO  03 
22  41 

0  1 

-20G0C 

ccoco 

D  i 

-  23  72  3 

-  23  7^" 

1 
i. 

1 

22  39 

22  ui 

01 
0  1 

-S3SS 

_0  CC  Q 

-0006 
-9391 

r  I 

n  i 

^r  23  7^  1 

to:  i" 

1 

22  '40 

0  1 

-  IS  33 3 

nn  it  l 

-9933 

D  ? 

S  9':  1 

gOCf 

6  3S 

C    70 

02 

999  3 

i  n  nn  5 

-  00  04 

G  ^ 

n  4 

g  gc  f 

£  37 
71  c2 

02 

04 

00  8 
00  0 

-9393 

-0009 

0  a 
n  a 

_p  pa  n 
_r  Cc  7 

71  5? 
7a  ar 

04 
0  4 

9397 

laoos 

00  00 

-93  9  2 

n  r 

10  p7  0 

74  41 
72  2^ 

04 
05 

0  00  9 
5993 

-9S93 
-0005 

r  r 

Ps5 

10  83  c 

ir  p7  'i 

72  25 
72  25 

0  5 
05 

20  00  0 
20  CC  7 

00  00 
-39  93 

10  83  7 
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^b 


05 


10  00  3       -3336 


66 


APPENDIX  II 

BLOCK  TRIANGULATION  PROGRAM 

The  Block  Triangulation  Program  performs  a  spatial  resection  and 
orientation  by  a  simultaneous  least  squares  adjustment  treating  photo  and 
survey  coordinates  and  the  elements  of  exterior  orientation  as  observed 
quantities.     As  significant  input  it  utilizes  values  of  the  camera  constant 
(focal  length);    corrected  photo  coordinates  output  from  COMORDCON; 
observed  values  of  exterior  orientation  equated  to  the  initial  approximations; 
and  various  job  keeping  information  such  as  number  of  photographs,  number 
of  unique  points,  photograph  number,  etc.     Additionally,  variance-covariance 
matrices  are  estimated  and  input  to  the  program  for  the  photo  coordinates, 
survey  coordinates,   and  the  elements  of  exterior  orientation. 

The  simplifying  assumption  that  there  is  no  correlation  between 
photo  coordinates,   survey  coordinates,  and  the  elements  of  exterior  orienta- 
tion is  employed.    It  is  further  assumed  that  no  significant  correlation  exists 
between  the  individual  survey  coordinates,  the  individual  elements  of  exterior 
orientation  and  individual  photographs. 

The  output  of  this  program  provides  the  minimum  variance  solution 
for  the  adjusted  values  of  the  elements  of  exterior  orientation  and  the  survey 
coordinates. 

This  solution  is  effected  through  the  following  mathematical  form- 
ulation: 

The  observation  equations  are  : 


V 

+ 

9        0 

B  A 

+ 

s     a 

BA 

+ 

c 

= 

0 

a 

V 

- 

e 

A 

+ 

0 

= 

0 

8 

V 

a 

A 

+ 

I 

( 

= 

0 
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In  matrix  form 


■ 

e      s" 

e" 

-      - 

V 

B  B 

A 

e 

e 

8 

• 

V 

+ 

-I    0 

Ad 

+ 

c 

s 

L.      J 

■ 

V 

0-1 

<r 

J 

- 

L.      J 

or 


V    +    BA+    €=0 

Then  by  exercising  the  method  of  LAGBANGE  for  realizing 
the  minimum  variance  solution  we  have: 

F   =   VTWV  -  2XT(V  +  B  A  +  O 


=    0 


thus 
(BTWB)  A   +    (BW?)    =    0 

where 


B   = 


EQUATION   A 


e 

B    . 
B 

e 

;      B  = 

^   R    ° 

^S.    a 

0         "Be.. 

8- 

B   = 

0        B3B 

8 

;bj  = 

ax, 


a  (Xq  ,  Y0 ,  Z0,  Xo»<Po>  k>o , 

ii 

,d(Xo,  Y0,  Z0,  Xo,<po»  ^o) 


bX 


d(x,Y,z; 

2^ 

La(x,Y,z; 


68 


w 


p                   -I 

w  w°- 
0        w 

- 

;      W    = 

0 

w2  ° 

• 

w3n 

;    Wj    = 

wxx    o 

0   Wyy 

\ 

0 

— 

a 

w  = 

Wi 
0 

a      0 

w2 

a 

;    Wj    = 

WyOyQ 

0                           Wyocoo 

~  a 

" 

8 

w  = 

0 

8      0 

w2 

\  8 

Wn 

;    Wj 

Wvv                   0 

o    Wyvwzz 

A  = 


a 

A 

8 

A 


BWe 


u 


a 

r-             "1 

6i 

a 

6! 

~6Xj" 

a 

A     = 

$3 

;       » 

8 

A    = 

52 

;      6;  = 

6Y, 
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• 

b 

66m 

e 

63 
_63a 

6Zj 

xi  -  xi 

"xr  -x," 

€       = 

o        «. 
Yi  -  Yi 

o          °° 
Yn    -   Yn 

» 

a 
€ 

00 
YI    -  Yl 

0          °° 
0        °° 

O               oo 

8 

0           CO 

Yi-  Yx 

<J            00 

Zi:  zl 

U                00 

Zan"    Z3n 

Ul 
U2 

» 

a 

Ul  =  BTwe 

8 

e    a 

-  We 

8     8 

U2    = 

=    B 

rwe 

-we 

with  m  =  the  number  of  photos  and  n  =  the  number  of  unique  points. 
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Expanding  EQUATION  A  provides; 


BT  -I     0 


BT     0  -I 


W 


W 


W 


B 

s 
B 

-I 

0 

0 

-I 

8 
A 

+ 

BT-I    0 

S 

s 

A. 

BT  0  -I 

and  after  multiplying  out, 


e            •             a 

BT  W  B    +    W 

•             s 

BT  W  B 

a              e 

BT  W  B 

8                     8                    8 

BT  W  B   +   W 

9 

A 

s 

A 


e  so 

BT  W  C  -  W  € 

8  as 

BT  W  e  -  W  ( 


0     EQUATION    B 


and  N  A   +    U    =    0 

9  9  9 

By  inspection  one  can  ascertain  that  BT  W  B  +  W    is  a  full  six  by  six  matrix  with 
one  photo,  but  for  m  photos,   a  block  diagonal  matrix  of  6  m  by  6  m  consisting  of 

9  8  8  8 

six  by  six  submatrices.    Analogously,    BT  W  B  and  BT  W  B  will  be  6  by  3  n  and 
3  n  by  6  respectively  for  one  photo  but  6  m  by  3  n  and  3  n  by  6  m  for  m  photos. 

8  8  8 

BTW  B  +  W  will  naturally  be  3  n  by  3  n  for  one  or  m  photos.     However, 

8  8  ■  I 

because  of  the  assumption  that  there  was  no  correlation  in  W  and  W,   BT  W  B  +  W 
is  a  block  diagonal  matrix  of  three  by  three  submatrices.     These  facts  concern- 
ing this  system  of  normal  equations  lend  EQUATION  B  to  a  simplified  method 
of  inversion  by  partitioning   [13]. 


N     = 


8                     8                    9 

l        e          B 

1 

BT  W  B   +   W 

BTW  B 

A        B 

8                    9 

8                    8                    8 

- 

__  _ 

BT  W  B 

!       BT  W  B    +   W 

C    |    D 

where  N  is  of  the  rth  order  and  A  and  D  are  of  order  6  m  and  3  n  respectively; 
r  =  6  m   +    3  n 


N" 


K     ]     L 


M 


N 


and       N"1   N   =    I 
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therefore 


It  can  then  be  verified  that; 


thus; 


N*  = 


A  K  +  BM  =    I 

A  L  +  B  N  =    0 

C  K  +  DM  =    0 

C  L  +  D  N  =1 

K    =  (A  -  B  D^C)"1 

L    =  KBD'1 

M  =  -D4C  K 

N   =  D^-D^C  L 


(A  -  B  D^C)* 
-D'C  K 


K  B  D"1 


D^-D^C  L 


e  a  a  f  t  •  c  » 

and  where     N   =    BT  W  B;         N   =    BT  W  B;         N   =    BT  W  B 


a  _ 


Kl  = 


C       8  S  C  -1 

1  mT 


c     a  s 


(N   +   W)  -  N(N  +  W)'1  N'  -K  N(N  +  W) 


8  8  C 

•(N  +  W)'1  NT  K 


j  (N  +  W)'1  -(N  +  W)"1   NT  L 


As  a  result  only  one  6  m  by  6  m  matrix,   and  n  three  by  three  matrices  must 
be  inverted  in  the  partitioned  matrix  to  effect  the  inversion  of  N,   an  rth  order 
square  matrix  (r  =  6m  +  3n),  for  m  photos  and  n  unique  points.   In  solving 
the  normal  equations  the  solution  then  to  the  alterations  to  current  approxi- 
mations is:    A   =    -N"1  U  and  the  adjusted  solution  vector  is  computed  as; 

Xa     =     X     +     A 

with  subsequent  iteration  converging  to  a  suitable  solution. 

Xa1  +  l=    X,,+    &t 

The  BLOCK  TRIANGULATION  PPOGBAM  is  capable  of  processing  up  to 
twelve  photographs,  36  unique  points,  and  a  selected  number  of  iterations 
with  but  minor  modifications.     Additionally,   it  provides  the  standard  error  of 
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unit  weight,   and  the  residuals,   standard  error,   and  the  variance-covariance 
matrices  of  adjusted  parameters. 

Other  than  standard  library  subroutines  this  program  utilizes 
SUBROUTINES  MATINV  and  COFET.     MATINV  is  merely  a  double  precision; 
matrix  inversion  routine  employing  a  bordering  technique;    COFEI  computes 
the  partial  derivatives  of  the  projective  equations  with  respect  to  the  elements 
of  exterior  orientation.    The  main  program  provides  for  updating  photograph 
and  point  information  for  use  on  successive  passes  of  the  iterative  process. 

The  terms  in  the  preceding  discussion  are  defined  as: 
V  =  residuals  on  observations 
B  =  partial  derivatives  of  the  projective  equations 
A,  5  =  alterations  to  current  estimate  of  the  unknowns 
c  -  functional  relationship  between  the  observations  and 

unknowns  evaluated  at  the  current  estimates 
X  =  Lagrange  multiplier 
W  =  weight  matrix 
x,y=  photo  coordinates 
X,Y,Z=  survey  coordinates 

-*o»    o»    o  >  _  eiements  of  exterior  orientation 

X,(p,OJ 


Superscripts 

refers  to  elements  of  exterior  orientation 

refers  to  survey  system 

refers  to  photo  system 

refers  to  observed  values 

refers  to  computed  values  using  approximations  to  unknowns 


none 

o 
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BLOCK  TRIANGULATION 
LIST  OF  PERTINENT  VARIABLES  AND  INPUT  DATA 


CARD  NO. 


4,5,6,7,8,9 


VARIABLE 

A 

Ul,  U2 

DE 

DS 

E 

R 

CC 

N 
WM 


IP 

Photo(l,l) 

Photo  (1,2-7) 
[X0,EE] 


10, 1L  12, 13, 14, 15  Photo(L14-19) 


16 

Photo(l,  1,1-2) 

17 

Photo(l,l,3) 

18,19,20 

Point(L  L  17-22) 

(X,  ES) 

21,22,23 

Pointy  L  8-10) 

PURPOSE 
Elements  of  N  and  N"1 
Constant  vector 

Alterations  to  elements  of  exterior  orientation 
Alterations  to  survey  coordinates 
Residuals  on  x,  y  photo  coordinates 
sin(x,  <p,oj)  cos(x,  (Aco) 

Camera  constant  (f  in  mm) 

Number  of  photographs 

A  scalar  multiplier  of  a  2  by  2  identity  matrix  (I) 
to  provide  an  estimated  variance-covariance  matrix 
for  observations  on  photo-point  coordinates 

Total  number  of  unique  points 

Number  of  points  in  photograph 

Current  approximation  and  observed  values 
(equated)  of  the  elements  of  exterior  orienta- 
tion (m.  and  radians) 

Estimated  variance-covariance  matrix  for 
exterior  orientation 

x,y  photo  coordinates  (COMCORDCON) 

Point  number 

Current  approximation  and  observed  values 
(equated)  of  survey  coordinates  (m) 

Estimated  variance-covariance  matrix  for 
survey  coordinates 


CARDS  3  through  15  repeat  by  m   number  of  photographs  and  16  through  23 
repeat  by  n     number  of  points  in  each  photograph. 

•    73 


MAIN 


"C  ""      "     bLUCK.    I Rl ANGULATION 
C UECK    LtUlbbblitlbbbbUUb 


0001  IKI'LILII     HLAL*f){  A-M.O-/I 

0002  ".  "  OlhFNSIC*    A1K10,  IB)  *  At2f  l»«  1*1  ,A22(15,15) 

0003  oim.'.suiK  iJUi(n»,n),A2/ni«;,ibJ,oe«UJ,osil5) 

0004  DIMENSION'   U  1  ( 1  a  1  •  112  {  1 5 )  t  Ah  f  18),PhUT0(3,4v),?ntNT(3,5,22) 
C005 UIMENSlON    XO  (  6  )  ,  C  £  ( t>  )  ,  hfc  (  6  ,  6  » 

"OOOo  "DlflCNS  I  0  n    OAT/.  (?)  ,  hlx,?  )  ,x'(  »)  , eS<  3)  ,«S«  3.  J) 

OOU/  OIKLNb  I  UN    L  (,*  )  ,  AllMr(fa,2>  ,U(2,-4)  ,K(6)  ,MONfH(  1  2  )  ,  I  F  1  1 1_  6  (  i  ) 

oooa  Dimjsia*  kli  iooi  ,kliuoo> 

0009  INJL'lilK    i>P,P,l  ,KEY(35)/35*0/ 

0010  '   OATA    MUNllH I ) / • J  AN. • / ,MON I H < 2 ) / ' FEB. • / ,MUNTh( 3 ) / 'MAR. « / .KONTH < 4 ) / • 
*AP«.,/t>0\rH(i)/«KAY     •/.MONTHlb)/*  JUNE  •  /  ,  HON  Jrlt  7  )  /  •  JUL  Y  •/  .MONTH  I  8  ) 

♦  /•aijo.  '/.MnNHiOJ/'SCPT'/fMONTMdOJ/'OCT.'/.MONTHdU/'NUV.'/tMONr 

*H<  m/'CLC.'/ 
OOU  "  DAIA    K.M,KN,KKJ/'N» ,  ».V  ,  •  t'/ 

0012  DAIA     IbLK/'     •/ 

001  3        ■"'     '  ~  1001  =0 

0014  DF=0.0 

0016  ICYcLE=0 

CO  I  7 IU  =  5 

0018  RfcAC(5.5lO)     HllLt.JObNUK 

0019  REAEI5, 2>>)    CCfNtlP 
C020 KrAU(5,29'>)     WM 

0021  ZTTT&rvi  ATI  F  10.3,?  T5) 

0022  UO    2C0    1  =  1, N 

0023  RcA0(5,603>    PhOTOC  I,  I) 

0024  P-HH(iru(  1  ,1  ) 

0025'  "     ULi    201    J  =  2,7  ": 

0026 KLAUli.500)    PIUHUI  I ,  J  )  ,  PHCHOC  I ,  J*6)  

0027  201  THGTOI  1  ,  J*6)  =P»iGTOC  I  ,J>'~ ~ 

0028  L2=13 

0021  "•                            LO   ?02    J  =  l,6 

0030  '    Ll=L2*l 

OOH  L2=L2*6 - " " 

00i2 202  KEAU(5,5C1>  I  PHUTUJ  I  ,K1  ,K  =  Ll  ,  L2)  

0033 Oil'  200  J=l,P  '••"" ' 

0034  KfcAUC>,jO?>     POINT  J  1,J,1),P01NI(  I, J, 2) 

0035  ftEAD(5,6G3>    P0!NT(I,J,3)  -  ■■  _     • 

0036  PUlNTt.  1  ,  J,  <.)-».»'. 

0037  " ~  POIivin  I,  J,  51  =0.0 

0038 POlNTt  l,J,6)=0.0 

0039  PGTNr<irJ,/r='KH 

0040  00  204  K=l7,lv 

0041  "'""■■  KEAfj(5,5U0)  PUI'4T(I,J,K),POINT(I,J,K»3) 

0042      204    PC1N1 ( I , J ,K*31 =PulNT ( i , J ,K> 

004 3  ~~~~^      ~~    """  L2  =  7 

_0044_- 00    205    KfJ.jrJ 

0u4  5  "Ll  =  l2+T  : 

0O46  L2=L2*3 


0047  205   RCAC(5,5C4)     {  POINT  nyj7U",L>Ll"YL2  T 

0048  200  CONTINUE 

_  . c . 

C LLEAK  NOKMALS ' 

"004'3  i<T=3*rP 

0050  Jl=b*N 

0051  4    ICYCLE=ICYCIE*I 

0052  00    5    1  =  1, IP 

0053  5-ii£YII)=0 

0U54 IOC    2    1*1, Jl 

"ooss  ui  n  y=o .  cr 

0056 00    I    JM.Jl 

C057  T  A11II,J)--0.0- 

0058  00    2    J=l,Kl 

"0059  2    A12(l,J)=0.0 

0060 UO    3    1  =  1, Kl 

"0061  IT2  T77VG .  C 

0062      OU    3    J  =  l,Kl 

0063  3"  A22(l ,J)=0.0" 
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C       "      CALL    PHOTO    DATA 
C 

-ooos  rpc=i 

C065  00    100    J-  l  ,  i4 

00<.6  LCC=13 

006  7  K  =  Ph.'JIO(  J,  I  ) 

OO08  ~~  "00    210     11=?, 7 

oooo  xf; c  i  i-i  )  =  ->rifir(!(  j,m 

00/0 "E £  (  I  F I  )=?nd I 1) ( J r~~ 

0071     00    210    JJ=l,6 

00/2    ~=  L0C*L0C*1 

00/3  210"  wF(  11  —  1  ,  JJ  )=^liLTU(J,LOC) 

00/4  »RITEU,601>    XO 

CO/5  *;<I  rt  lo.o'/l  )    HE 

"  CO  76 50T~FfUVAT  ( 6£1 5 . 5  ,  /  (6X157711  ' 

00/7  K(  I  )=USl.\(XO(<.  )> 

0073     7  R(2)*D„IMX0(i>>)  '"    '* 

00/9  K(  ))=uilii(XO(o)  ) 

OObO  k<4)=DC0SUCU4>> 

00b  I K(t>)=uCUS(XU<5)  > 

_oob2  s"(6T=ijcai(xo(6n 

c 

"»i C ADOING  TERMS  JOE  TO  CONSTRAINTS  QW  EXTERIOR  ORIENTATION  tLEXtNTS 

C_ 

0083"  L=o*(J-l)  """' 

00b4 UO  10  11  =  1,6 

"C0b5  IFUrl  I  1.1  I  r."N_"."0.~OJ  OF=uF*1.0 

0086    VUF  =  VUK»V,E(llf  IM*(CE{IM-X01II))**2 

0087  DO  10  JJ=l,6  

OObb  U1(L*1 I )=0l (L*1I )-Wt(Il,JJ)*IEC(JJ>-XO(JJ)» 

0089  10  All(L*I l,L*JJ)=All(L*Il,L*JJ)*tf£(ll,JJ)        ~ " 

C 

C      CALL  POINT  DATS 

C 

0090  "       00  100  1=1, K 

0091  l)AU(  I  )=r\irjl  (J,  I, I) 

0092  """',    0AlA(2)=PUlNr(J,l,2>  "  

0093  P=PCINT( J, I, 3) 


0094  LGC  =  3   "  "~"  

0095  UO  215  11  =  1,2 

0096  '    _0  215  jJ=l,2 
C097  LOC=LCC»l 

0096  215  wm,JJ)  =  POIUT{J,  l.LOO 

C099  U(J  21b  11  =  1,3 

oioo  on  ?\u   jj=i,3 

0101  LIK  =  L0C*1 

0102  — —  —  216  *S<  Il.JJl-POlUTUtl.LOC) 

0103  00  217  1 1=17, 19 

0104  X(II-l_)=PfJHTU,If  II) 
0105 217  liSt  H-16).=PIUNT<J,I,  11*3) 

0106  nRlTf (o,o01)" X 

0107  wRl TH(6,o01 »  ES 

0108      —   0F  =  0F  +  2.0 

0109  LL=3*(P-1> 

.  .  C  "" 

C      CALL  PARIIALS 


0110  CALL  CCFt I (X,K,0,XO,CC,XC,VCI 

.,  ; c 

C  •      COCPUIF    NORMALS 

0111  ITcSt  =  l  —■    ■   -  — 

0112 IF(KtYlP).GT.O)    GO    TO    12 

0113  Kfc Y ( P )  =  P  

C114  IT£ST=0 

0115  12   00    15    II=t,6  '"" 

0116  00    15    JJ  =  l,2 

0117  ATEMi»(  I  I,JJ)=0.0   ""•"    '"  ■"""" 

Ollb 00    15    KK=1,2* 

0119  IT~AT£Mf  (  I J  ,'J  J  >  ="A TcHPTl  I",  JJ  f+O'liCXTIT  J  **  ( TfXCJST 

0120  L (  I  )=UArA( ll-XC 

0121  L(2>=0ATA(2)-YC 

0122  </0F=V0F**(  l,l)*E(l)**2*W(2,2>*E(2)**2 
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UU'3 
0124 
0125 
01  2ft 
0127 
0128 


coMfuri?  ah 
on   W   1 1- 1  ,<» 

Oil    lo    KK=l  ,2 

tl  (L*l  I  >=Ul  ( L  ♦  I  l)»ATb.*P<  I  l,K.O*£(KK) 

0(1     I  /    J.J-  I  ,6 

00     17    KK=1,2 


16 

17All(l«Il.L*JJ)=AllU+ll,LOJ)*ATEMP(lI,KiO*aUK,JJ) 


0129 
0130 

01  31 

"o  inr 


~Yf 


COMPUTE    A12 

00    19    I  1  =  1,6 

DO    \'>    JJ=7,9  ""•" 

JO    19    isK  =  l,2_    _        _  _. 

"  A  I  2  f  L"  ♦  1  1  ,  LL  ♦  Jj-6  )  "AT2"  ( L* 1 1  *\.L+S3='bT*iFC&KV  (II. K'kY*&  fit* ',  J JT 


CUKPUfC    A21 


"01  33" 
0134 
013r> 
0146 
0U7 
CI  3d 

"0139" 
0140 
0141 
0142 
0143 
0144 


■  23 


"26 


CUrtPurK    A22  

tjc  ~i  3   i  r=  7  ,  9 

uO    23    JJ-1,2 

AIci'tM  1  1-6, JJ)=0.0 

00    2  3    KK=l,2 

A  TEMPI  1  1-6,  JJ)=AICMJ>{  I  1-6,  JJJ+QjKK,  II  )*W(KK,  J  J)     ' 

00    21    11=1,3 

h»LL>iT 

00    2o    KK=l,2 

U2(M}=02(M)»AlEMPUl,KK)*E<KK) 

OU    2o    JJ-7,9 

A22(f-,LL»JJ-6)=A22{M,LL  +  JJ-6)»ATEMiMII,KK.)*BUK.,JJ) 

in  m-sr.EO.n   go  to  27 


0145 
0146 
0147 

"ClVa" 

0149 
0150 
0151 
0152 


adding  it-u'.s  oue   ro  constraint  on  survey  coordinates 
rfOFsvof**s(  ii,  ii  )»(ts(  in-xmi  >»»2 

1F(»S(  1  1,11). NE. 0.0)    OF=OF*1.0 
00    2M    JJ=l,3 


A  22  "f  N  ,  L  L  ♦  J  J  YZATnH  ,  I L  *JJ)*WS"(  II, Jj) 
2d    U2(X)=U2(MI-wS(  Il.JJ)*(ES(JJ)-XUJ)) 
2  7    COM  I. SUE 
100   CONTINUE 

"1F(  lUUI.cO.l)    CO    IC300 


~01>3 
0154 
0  Ij'j 
0156 

"or;  7" 

015H 
0159 
0160 
0161 
0162 
"01 6V" 


C 
C 

C 


0164 

0165 

0166 

0167 

oioa 

Olo9 

0170 

01/1 

0172 

0173 

FORKING    Oil 

Call-  ma r l nv I A22i Ki) 

00   IS    I  *  1 , K I 

00     »•>    J»l,Jl 

_A22TI  I  .JJ^O.O 

DO    R   K«] ,Ki 

35  A22TU,J)=A22TU,J)»A22U,K)»A12{J,K) 
00    36    1  =  1, J I  "       " 

00     36    J=l,Jl 
00     36    iv'lfKl 

36  AU<I,J)  =  Alt(t,J) -A12(I,K)»A22T(K,J) 

Call  Hat  RV'f ATI73T1 

FORMING  612 

"   00  37  1  =  1,  Jl 

00  37  J=l,Kl        

ol2f{ 1 ,J)=0.0 
DO  3/  K=l,Kl 
il    lil2lll.JJM>l2l(I,J)*Al2<l,K)*A22«K,JJ 
OU  3«  1  =  1  . J  1 
00     in    J  =  l,M 
A12( I , J) =0.0 


3a 


uo  3-j  n  =  i;jt. 

Al2ll,J)=Al2(I,J)-AlllI,X)*612t«K,J> 
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FORCING  H22 


0174 

01/6 
017/ 


01  /tt 
01/9 
OlHO 
0181 
01t2 
01BJ 
01&4" 
0185 
0U6 

oia/ 

GU8 
0189 


GO 

39    1=1, Kl 

00 

39    J  ="l  ,  K  I 

00 

39    K=l,Jl 

39 


A22(l,J)-A22(I,J)-A22TU,K>*Al2(K,J) 
SOLVING  NOrtKALS 


LIU    40     1  =  1  ,  Jl 
OE(l)=0.0 
OlJ    ',1     J=  I  ,  J  1 

41  Otll)=0tlll-AIK1,JMUI|J) 
00    tO    J  =  l,M 

JHQ_L'FI  1 )  =UE  (  1  1-A121  I,J)*U2(J) 
U.J    42     l=l,Kl 
0S(  1)=0.U 
60    43    J=l,Kl 
A3   GSt  I)=OS<  I  I-A22I  I,J)*U2U) 
00    42    J=l,Jl 

42  US(  l)=0Sll  )-A12(J,I)*UHJ) 


0190 
01  yl 
CI  1/ 
0193 
0194 
0195 
C19o 
019/ 

019b" 

0199 

0200 

0201 

0202 

0203 

020* 

0205 

020b 

020/ 

020(3 

0209 

"0210" 
021  1 
0212 
0213 
0214 
0215 

10216" 
0217 
0218 
0219 
0220 
0221 

"0222" 
0223 
0224 
0225 
0226 
0227 


101 


115 


UU 
K  =  k 
01) 

CUiN 
U  (j 
j>  =  p 

uo 

>tJl 

Co.* 

in 

DF« 
VOI- 

^_JUG 

00 
3C0    CAL 
C.AL 
UF  = 
VArt 
UNS 
\.AL 
00 
UU 
All 
JO 
0(1 
"A  2  2 
NOF 
00 
IF( 
l  r  i 
uo 

"IF  I 
1  II 
I  I'M 
1  =  1 

J  =  J 

MR  I 


306 


3u7 


301 
302 


iOi 
304 


0228 
02^9 
0230 
0231" 
0232 
"02  i  3 
0234 
02  35' 
02  36 


»F  ,0 
wRl 

rt  ►!  ( 
/Kl 

"CO 

16  = 


uU 
X01 
"30i»    cLl 
WKI 


LY 

115 

lull 
101 
TOl 

nw 

Ul  1 
115 

;n( 

TIN 
ICY 
0.0 
=  C. 
I*! 
Ill " 

L  M 
L  M 
If— 
=  V0 
TOR 
L     1 

l'J6 
300 
(I. 
30/ 
30/ 
(I, 
«0f 
301 

I  1  1 
(•£  = 
303 

II  1 
l'.t  = 
L  =  F 
-I 
-t 
TFJ 
Nil 

IC( 

It  ( 

Icl 
3iO 
o«( 
»05 
J)  = 
J)  = 
Td 


t-U 


ALTERATIONS 

1  =  1  ,N 
'J  I  I  ,  1  > 

l/  =  2,  / 

i ,  m=.»Horo(  I,lZ)*Ud6*(  I-UMZ- 
LF 

J=l,K 
I  I  1, J. 3) 

U=17,19 

I  ,JYl2»  =  k»0lnn  I  ,J,  IZ  >  ♦OSC  3»~(>-l  )»U- 16 f 

OF 

CLC.LI.6)  GOTO  4" 


AT1NV1A11 ,Ji J _ 

A  IIMVIA22.K1) 

Jl-Kl 

F/liF 

cOSQRUVAR) 

IJAI  1  "1  (  I  YEA*',  IKONI H.YuAY  ,  1 1 IME  f 

1=1. Jl 

J=l,Jl 
J)=A1 1  I  I , J)*VAR 

1  =  1  ,Kl 

J=l,Kl  '    • 


J)=A22( I ,J)»VAR 


1=1,24 
Mc.Ll.  16C000)    OOJTO    302 

lTll,t-3oC000 

J =1,59        

PF.Lf  .6"000">~G~U  TU  304 

I  IIKF-oOCO 

LO.MI  I  I  l^tl/100.0 


WW 


6.J000)     11  I iLt, JOBNUM,  IOAY.MGMIhl I  MONTH), I  YEAR , I , J , TIKE ,N.NO_ 

OK 

1,31 

2,  iC 

o,  J( 

1  = 
i-i; 

j  =  i 
USv 

t'HUl 
6,  i(. 


*JOi) 
1005) 

>001)_  

<I,N 

) . 

1,6 
„<U(A1  It  16*J,  I6  +  J)  ) 
;rc(l  ,J*7)-fH0I0tI,J*U"  "" 
J002)     l.l-'MCTCl l,J»l),J»l,6>,X0,EE,  t»»HOTO(I,J> , J" 14,49,71 
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0217  t.R  I  If  I  /,,   1007)      (  (  Al  1  (  I6»K,  [6»J  )  ,  J«  1  .6)  ,K»  1  ,6) 

02  til  Hill  It  Ih,  IUOVI 

0239  R(l)-OSIN(PHuro( 1,5) ) 

0240  K(/)M)SIN(PnfiTU(  1,6)1 

0241  R(3)=0SlN(PHCTO( I,7>1 

0242  R(4)=0C0S(PHUT0<  1,5)  ) 
024J  K(5)=DC0S(PH0T0(  1.6)  ) 
024<,  Klol-iiUJjIPnutull.ni 

0245  K=pnuroii,n 

u24o  uo  i.">  j-i,k 

024/  P=P01NT(I,J,3) 

02<>8  KK=3*(P-1) 

02^9  ox  =  puini  ( i ,  j,  m-PHuro(  i,2) 

0250  UY  =  t'UlNI  I  I  ,J,  18)-HHUT()(  I  ,  1) 

0251  02  =  P01NT(  I  ,J, 19) -PHOTO!  1*41 

0252  XT  =  uA*i<<5)*R(4)*DY*lR(6)»ft(  l)«R(3)*R(2)*R(4))*0ZMft(3)*R(l)-R(6)«R 
*<2)«R(4) ) 

0251  YT=-UA*R(t>)*R(  1)»0Y*(R(6)*P<4)-R(3)*R(2)*R(1))*D:»(R(3)*R(4)»RC6)* 

«R(2)*R(1) I 

0254  ZT=DX*R(2)-DY»R(3)*R15)*DZ*R(6)*R(5) 

0255  XC=CC*XT/ZT 

0256  YC=CC*YT/2T 

025/  EU)  =  P0INT1  I,J,l)-XC 

0258  E(2)=PCINT< I,J,2)-YC 

0259  «RI  IEU.3004)  t.P.POINTCI.J, t ),POINT( I,J,2),E(l),E(2) 

0260  ,         IFUEY(P).EO.O)  GO  10  325 
G2&1  K£Y(P)=0 

0262  00  320  JJ=l,3 

0263  X(Jj)=DS'jRT(A22(KK»JJ,KK*JJ)) 

0204  320    fcS(  JJ)=POINT<  I  ,  J  ,  J J  ♦  19 ) -PU I  NT (  1  ,J,  JJU6) 

0265  wRIIE(2,3006)     P,(PUINTII,J,JJ),JJ=17,19),X,ES,(PU1NT(1  ,J,JJ),JJ  =  8. 
•16,4) 

0266  URI  IE ■  (2,300  R)  I  IA221KKH  1  ,KK*  J  J  )  ,  J  J=  I  ,  3)  ,  1  I*  t  •  3  1 
026/                WRITEI2,3009) 

0268  1F( J/2*2.EQ.J)  Wtt I TE ( 2, 3010) 

0269  325  CONTINUE 
02/0  330  CONTINUE 

02/1  3000  FORMATl  •  1  '  ,26(/)  ,T<.7,  7A<.,/T53,  •  JOB  NUMBER ' , 16, /T53, 'DATE ' , I  3, 1 X,  A4 

♦  ,15,/154,'UME', 13,' :',I2, •:• ,F4.1,/f51, 'NUMBER  UF  PHOTOS  «  ',12,/ 
*149, 'DECREES  UF  FREEOGM  = • , I  5 , / 144 , • UNI  T  S1ANDAR0  ERROR  *  •,012.5) 

0272  3001  FORMAT (• I • , T58 , • Rt SUL I S ' , / 1 5 1 , • EXTERIOR  ORIENTATION') 

0273  3002  FORMAT!//, T10, 'PHOTO  NO.  ',I2,T27,'Xn  ( ME TERS) • , T41 , • YO  (METERS)', 

*T57,'ZG  (METERS) ', I  73, 'KAPPA  UAO. ) ' , 191 , • PHI  ( RAD. )', T 107 , 'OMEGA 

♦  (RAO. ) '//I20,3F16.3,3D17.6,//T9,'ST0.CRR0R' , T 20, 40 1 6. 4 ,2 ( IX, 016. 4) 
*,//T9,'RESI0UALS',I20,40l6.4,2(lX,Dl6.4),//T9,'WEIGHTS',T20,6F16.3 

♦  ,//) 

02/4  3003  F0RMATIT58, 'RESULTS' ,/T53, 'PHOTO  C00RD1NA TES ' , /T4Z, • ( ALL  WEIGHTS  T 

♦AKEN  AS  ',F7.l,' )',//T27, 'PHOTO  NO. ', T4i , 'POINT  N0.',T55,'X  (MM)*, 

02/4  )0U'.   (  ilMAI  (  'd',  Ill),i?,T'.4,|  l,(')0,?ni.  1,  T7J,  ?D1 4.4) 

_02  75 3J05  FOKKAr  I  T  ,",  '.U.SULIS'  ,  /  152, • SUKVLY  COdROlNAl ES  •  ) 

C2/o  ,  3uGo    f LMIAl  (//( 31 . 'POlNl    NO.     ' , I  3 , Tv2 , • X • , T68 , • Y » , T84, • Z ' , //T42 . 3F l&.Vt 

*//!  J/,  •',  ID.    IKKiJR' ,T4H,li)l6.4,//T32,  'RCS10UALS'  ,  (  43,  3D16.4.//T32*  ' 

*.,i   I  (,lil  •,'./.»!   lo.  »,//  I 
027/  300/    FURFAI  (  'u'  ,  I4H,  WARI  ANCE/LOV AR I ANCE    KATR  I  X  •  //  (  T  16,  6E  I  5  .  5  ,  /  )  ) 

0?/H  "  3u(lli    rii.'.IM  I  'D',l',a,'VAl<UNC(7CnVARUNCC    MAIRIX'//(T36,3£IS.5»/)) 

02/9 3e09    f  OKf-'Al  ('()',//) 

02nO  3ulo' KikHAl  I  •  I  '  ,7) 

0/t>l  >0(l    FlIRf  Al  (?F  10. 51 

02o2  501     FUKMAl  (6H0.5) 

0283  502"  FilRMU  (2F10.4) 

0284  :>03   FCiKllAl  1  15) 
028j  s04    FOKFAI (  IF  10.3) 


O/o/,  •jO'j  ►"OCMAM  '  iWiUTtl    NUMB t R  •  ,  I  3  ) 

028/ 506  FUXMAI  (  'OPQI.Sf     .IUMBcR  '  ,  14, /_)_ 

"028S    "  510  FURMA1 ( ZA4.2X, 14) " 

02o9 600  FORMAT  I IX,D15.7)  

0290"  60  3  FGP.MATIFi.O)-    

0291 STG*        

0292  DEBUG"- SUbCHK 


0293  END 
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MAUNV 

00</l  Mli'i'i.ur  I  <T    MAT  |.-4V(  A(NI 

C002  IMl'Lltll     Kt'AL*H(  A-H,0-Z  1 

0001  D  I  KHNSTO  r  A  (N,  N  J'.Ti  (  ftO  JVC  (1*01* 

COO'.  M  =  N-1 

C005  4(1,11=1.0/4(1,1)      

COOb  lF(l')     ?9o0,2r/0o,290O 

000/  2<V0G    uU    i<yO'j     1  =  1, is 

COOK  L  =  I ♦ t 


GOO'/  '           DO   2001    J-V,  I 

0010  UJ  1=0.0 

0011  "  "        . 2901    C(J1=0.0 
001/  i,<)    ^  ;<J2    Jrl  ,| 
C013  "■••  ■""           "'  00    2-702    *.*l,l 

00 1 <, i">  U  1  =  I'  (  J  I  -  A  (  K  ,  J  )  ♦  A  ( l ,  K) 

CC 1 5  2"502~  C"(  iU  »C  f  K 1  -A  ( K ,  J  V*  A  TJ"i  LI 

0016  0*A(LfLl 

0017  ~  ~'-     DO    200.1    J=l,l 
001H  2903    C=D«C< J1*A(L,J) 

0019  U=l.O/U 

0020        CO    2?0<>    J=1,I 

0021  MJ,L)=CU)*0 

0022  AU,J)=B(J)»0 

0023  0U'270A    K=1,I 

CG<:'.  2vC'a  (j^)=A(j,K)*b(K)*ClJ)*0 

0025  2-JG5   A(L.Ll-U 

C026 2<i06   KEtuKN  

'0027  Ztfb 
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cum 

0001  iUoauUTI-**-    CuFfcHUfRthiXO.CC. XC.YCJ 
C 

c  '  '    coci'urcs  i>  n><  L<rtKini<  ano  intckio*  liekcnts'inciuoing  cc' 

C  HHUl^rl)    '.mOHK     (  X,Y, Z,K,P,W) 
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iy  f*  (oy*cw»cp*u/*sw*c»>)  > 

"0033  YC=CC*YT*( 1.0/ZT)  
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